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ABSTRACT 

We report the discovery of significant localized structures in the projected two-dimensional (2D) 
spatial distributions of the Globular Cluster (GC) systems of the ten brightest galaxies in the Virgo 
Gluster. We use catalogs of GGs extracted from the HST AGS Virgo Cluster Survey (ACSVCS) 
imaging data, complemented, when available, by additional archival ACS data. These structures 
have projected sizes ranging from ~ 5" to few arc-minutes I to ~ 25 kpc). Their morphologies 
range from localized, circular, to coherent, complex shapes resembling arcs and streams. The largest 
structures are preferentially aligned with the major axis of the host galaxy. A few relatively smaller 
structures follow the minor axis. Differences in the shape and significance of the GC structures can be 
noticed by investigating the spatial distribution of GCs grouped by color and luminosity. The largest 
coherent GC structures are located in low-density regions within the Virgo cluster. This trend is more 
evident in the red GC population, believed to form in mergers involving late-type galaxies. We suggest 
that GC over-densities may be driven by either accretion of satellite galaxies, major dissipationless 
mergers or wet dissipation mergers. We discuss caveats to these scenarios, and estimate the masses 
of the potential progenitors galaxies. These masses range in the interval 10® ® —larger than 
those of the Local Group dwarf galaxies. 


1. INTRODUCTION 

The large scale distribution of GGs has been studied 
in depth in several elliptical galaxies. Radially, GCs are 
known to extend to larger radii t han the diffuse stellar 


light of the galaxies (Harris 1986), with the most lumi¬ 
nous galaxies having rela tively more centrally c oncen- 


trated GCs distributions ([Ashman fc Zepf 1998). The 


radial profiles of red (metal-rich) and blue (metal-poor) 
GCs differ, with red GCs matching more closely the light 
profile of the stellar spheroid - except for central flat¬ 
tened cores in the GC distribution - and blue GCs pre- 


seating a flatter extended distribution (Rhode & Zepf 

20011 IDirsch et al.||2003| [Bassino et al. |20U6t 

iVIineo et 

al. 2014). Studies of the 2D distributions of 

GC sys- 


consistent, on average, with those of the stellar spheroid 
of th e host galax y (e.g. NGC4471 in Rhode fc Zepf 
[2001] NGC1399 in|Dirsch et al.|2003||Hassino et ai.|2UU6l 
NUC 3379, NGC44 06~lSd NGC4594 in jkhode fc Zip 
2004| NGC4636 in [Dirsch et al.||2005| multiple galax- 
iHmjrs fc Rhode||20121 NGC 3585 and NGC5812 


les m 


Lane et al. 201^ NGC4621 in [Bonfini et al. 2012 


and D’Abrusco et al.||2013). However, more recent work 


using the large sample of ACSVCS galaxies (Wang et al. 


2013) has uncovered highly anisotropic azimuthal distri¬ 


butions of both red and blue GCs; these GCs are prefer¬ 
entially aligned with the major axis in host galaxies with 
significant elongation and intermediate to high luminos¬ 
ity. 

These characteristics of the GC systems provide in- 
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sight on the evolution of the host galaxy, as well as on 
the formation of GCs. The relative importance of dissi¬ 
pative merging or of general galaxy collapse has been 
debated in the literature, mostly based on analytical 


models and observational inferences ( Searle fc Zinn|19T8 
iman&Zegf 1992 and references therein; [Forbes 


see 


As 


eFaT 

smooth 


1997 Bekki & Shioya 1997). Recently, N-body 


red pc^cle hydrodynam ics simulations (Bekki fc] 


Shioya|19M Bekki et al.|2002| ) have reproduced the gen¬ 
eral radial spatial distributions and metallicity properties 


spatial distributions and metallicity properties 
of red and blue GC systems (but not the observed ratio of 
red and blue GCs) under the assumption that the ellip¬ 
tical galaxy was the result of major dissipative merging, 
as long as the merging event happened early enough in 
the life of the galaxy. Tidal stripping and GC accretion 
in interacting elliptical gal axies could help rec oncile the 
simulations with the data (Bekki et al. 2003). Dissipa¬ 
tionless merging of spirals could also give rise to th ese GC 
systems, but should retain kinematic signatures ( Bekki] 

etar][2005t . 


while the results discussed above still do not provide 
a unique scenario of GC and galaxy formation, it is clear 
that both major merging (dissipative and more likely dis¬ 
sipationless) and tidal interactions play important roles. 
More recently, localized 2D inhomogeneities and streams 
have been found in GC systems of elliptical galaxies, 
emphasizing the role of minor merging and companion 
galaxy tidal disruption and accretion. These results in¬ 
clude a spatially extended spiral-arm shaped over-densiW 


et al. 

2012 

D’Abrusco et al. 

2013D, 

anisotropies in l\(J(J4t)4y anc 

rTTUC 

al. 

2014ajb 

. These features 

were c 


and more localized 


using uniform and complete GC samples from Hubble 
Space Telescope (HST) observations. At larger galac- 
tocentric radii, streams in the GC distributions have 
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been observed with gro und-based c overage in M8 7 (Ro- 
manowskv et al.||2012[), NGC4365 (iBlorn et al.||2(jl2b| 

and NGC4278 ( |D^XE?bsco et al.||2Ul4b|) - 

In this fourth paper in our series ( [D’Abrusco et ah 


photo metric properties and spatial extent (Jordan et al. 


2013 2014a|b|, we investigate how common these 2U 


anisotropies and streams of GC may be, by studying 
the 2D GC distributions of the ten brightest galaxies 
in the Virgo cluster observed in the ACSVCS survey 
(Section]^. This sample provides a unique catalog of 
GCs in tne central regions of massive Virgo early-type 
galaxies. As in our previous works, we analyze the 2D 
distribution of GC systems in elliptical galaxies with a 
K-Nearest-Neighbor (KNN) technique supplemented by 
Monte Carlo simulations. 

The Virgo cluster is close enough to permit the de¬ 
tection of GCs with HST imaging, down to relatively 
low flux levels with limited integration times. While the 
dense central environment of Virgo may discourage ac¬ 
cretion of companions, Virgo provides a range of local 
densities for its galaxies. Even in its densest environment 
there is observational evidence of accretion of a compan¬ 
ion galaxy, su ggested by the kinematic ally detected GC 
shell in M87 (Romanowsky et al.||2012 ). With our sam¬ 


ple, we can therefore begin to investigate observationally 
the past accretion history of galaxies in a range of envi¬ 
ronments, using the GC population as a potential marker 


of these events. Moreover, our method (D’Abrusco et al. 


20131 is more accurate for larger samples of GCs, such 


Virgo cluster ((Harris & van den Bergh 

1981 

IBrodie & 

Stradei'l 2UU6 Peng et ai.||2l)U8 Ceorgiev et al. 

2010 I). 


Section [3] we describe the data and the method. In Sec¬ 
tion]^ we present the results of the analysis of the spatial 
distribution of GCs of each individual Virgo galaxy. In 
Section we discuss the results obtained, focusing on 
the spatial distribution of the GCs structures, and we 
infer the properties of accreted companion galaxies that 
could have originated the observed GCs over-densities. 
We summarize our results in Section El We have used 


galaxy distances from Blakeslee et al. (2009). 


2. THE DATA 

We used data from the observations of the ten bright¬ 
est galaxies in the HST Advanced Ca mera for Survey s 
(ACS) Virgo Cluster Survey (ACSVCS|C6te et al.j2004|), 
which observed the central regions of the 100 brightest 
early-type galaxies in the Virgo cluster, although with 
incomplete spatial coverage. For two galaxies (NGC4649 
and NGC4365) we were able to use additional HST data 
to cover the full area of the galaxy. 

The ACSVCS galaxies were observed in two bands, 
with the F475W and F850LP HST ACS filters, corre¬ 
sponding to the g and z Sloan Digital Sky Survey filters. 
Seven of the galaxies in our sample were covered with 
a single central ACS pointing; the NGC4472 observa¬ 
tions were arranged in a triangle-shaped three-pointing 
mosaic reaching the D 25 elliptical isophote along the ma¬ 
jor axis. NGC4649 and NGC4365 were covered with a 
single pointing in the ACSVCS, but subsequent observa¬ 


tions produced ful l ACS coverage (Strader et al. 2012_ 
Blom et al.||20 12b). The GC cata logs tor our work were 
extracted trom|Jordan et aLl(|2009|), which includes 12763 
bona fide GC candidates, iaentihed on the basis of their 


2004). 

The number of GCs in the ACSVCS catalog associated 
to the ten brightest Virgo galaxies is 705 3, representing 
^ 55% of the total | Jordan et al.| (|2009|) catalog. This 
number represents the traction of the CC population of 
these galaxies that could be detected in ACSVCS obser¬ 
vations. The basic information about the ten brightest 
galaxies in the ACSVCS list of targets are summarized in 
Table 1. Visual magnitudes Bt, numbers of GC candi¬ 
dates an d optical morphology are all based on ACSVCS 
data (see|C6te et al.|[2004 Jordan et al.||2009 Ferrarese 


|et al.||20L)6[ tor Bt magnitudes, CC numbers and mor¬ 

phology, respectively). The details of the additional data 
used for NGC4649 and NGC4365 are given below. 

NGC4649 : a joint Chandra-BST program (P.I. Fab- 
biano) completed the single central ACSVCS pointing 
with a mosaic of five additional ACS pointings in the 
same F475W and F850LP filters, with comparable inte¬ 
gration times. The final mosaic covers the whole area 
of the galaxy comprised within the D 25 (~4.5 re), with 
the exclusion of the region occupied by the companion 
spiral galaxy NGC4647 (i.e. excl uding ~ 15% of the 


NGC4649 area). From these data, Strader et al. (2012) 
extracted a catalog of 1516 GC candidates that was used 
to charact erize the spatial distrib ution of the GCs in 
NGC4649 (D’Abrusco et al. 2014a). 


NGC4365 r Blom et al. 12Dl2b) used six separate ACS 
pointings in the same F475W and F850LP (P.I. Sivakoff) 
publicly available on the Hubble Legacy Archive to com¬ 
plement the single central ACSVCS pointing, and pro¬ 
duced a catalog of 2163 GCs covering the entire galaxy 
out to ^ 1.3 D 25 ("^3 Te). The magnitudes an d sizes 
of the GCs in the central pointing measured by |Blom| 


(|2012b|) are compat ible with the same properties 
Jordan et al .| (|2009|). Here we use the catalog of 
BIc ' 


et al. 
from 

GCs extracted by |Blom et al. (|2012b|). 

For each galaxy in 'Fable 1, we investigate the spatial 
distribution of the entire sample of GCs for each galaxy 
and of color and lu minosity classes, f ollowin g the pro¬ 
cedure outlined by D’Abrusco et al. (2013). The red 
and blue GC classes were determined by htting the his¬ 
togram of the g—z color distribution with a double Gaus¬ 
sian model, and choosing the threshold corresponding to 
the 50% probability for a GC to belong to the red or 
blue class. Although for some of the galaxies the color 
histograms are better fit by a single Gaussian model, 
we adopted the color threshold based on the double- 
Gaussian ht to define the color classes in all cases. 

For our analysis, we have also considered a split of the 
GC samples in luminosity using, for all galaxies, the fixed 
magnitude threshold g < 23 mag, corresponding to an 
average absolute magnitude < ranging between 

^ —7.9 and ~ —8.3 Mag (see Table for the absolute 
magnitude threshold for each galaxyp These luminos¬ 
ity thresholds produce high-L and low-L GC subsam¬ 
ples of similar size for all the galaxies in our sample. 
While setting the brightness thresholds as the turnoff lu¬ 
minosity on the GC luminosity function may be more 
physically meaningful, as it would correspond to a mass 
threshold, that choice would produce highly unbalanced 
high-L and low-L classes, making our analysis statisti¬ 
cally unreliable. The definitions used to define color and 
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TABLE 1 

Properties of the ten brightest Virgo cluster galaxies. 


Galaxy 

Bt 

(a) 

Dl [Mpe] 

(b) 

% D 25 area 
covered 

(c) 

iVxot 

(d) 

[N.ed.lVblue] 

[Nu 

ghL) .Nlowl] 

(f) 

Optical 

morphology 

(g) 

r. ["] 

(1.) 

D 2 t [re] 

(i) 

A(R.A.) 

a 

A(Doc.) 

S) 

w 

j\^nresQ Mlin,(GCs 
[Mag] [Mag] 

(□) (P) 

NGC4472 (M49) 

9.31 

16.2 


1206 


708,498 

[751,455] 

E2/S0'l-(2) 

209.15 

2.9 

- 0.004°(14.7") 

- 0.004°(13.4") 

1.14 

-8.10 

-6.4, -7.3 

NGC4486 (M87) 

9.58 

16.8 

-40% 

1745 

[1006, 739] 

[698.1047] 

EO 

171.71 

2.9 

- 0.002°(6.1") 

- 0.002°{7.2") 

1.21 

-8.11 

-6.4, -7.3 

NGC4649 (M60)^ 

9.81 

16.6 

~40%(~85%) 

807(1603) 


841,762 


[1136,467] 

S0-l_(2) 

99.34 

4.5 

- 0.005°(18.6") 

- 0.004®(15.7") 

1.18 

-8.09 

-6.3, -7.2 

NGC4406 (M86) 

10.06 

16.4 

-40% 

367 


140.227 



218,149 


S0_l_(3)/E3 

411.84 

1.3 

- 0.005°(18.6") 

- 0.004°(15.7") 

1.15 

-8.26 

-6.5, -7.4 

NGC4382 (M85) 

10.09 

17.1 

-45% 

507 


283,224 



241,266 


SO 1 (3) pec 

170.82 

2.5 

- 0.002°(6.7") 

- 0.002°{5.8") 

1.14 

-8.26 

-6.5, -7.4 

NGC4374 (M84) 

10.26 

17.0 

-45% 

506 


251,255 



296,210 


El 

149.73 

2.6 

- 0.002°(7.1") 

- 0.002°{6.4") 

1.12 

-8.36 

-6.6. -7.5 

NGC4526 

10.61 

15.3 

-75% 

244 


134,110 



124,120 


S0-3-{6) 



- 0.002®(6.9") 

- 0.002°(6.1") 

1.10 

-8.14 

-6.4, -7.3 


10.51 

22.0 

-45%(100%) 

907(2163) 

[1149,1014] 

[883,1280] 

E3 

116.14 

3.6 

- 0.005°(17.6") 

- 0.005®(18.9") 

1.16 

-8.82 

-7.1, -8.0 

NGC4621 (M59) 

10.76 

15.7 

-75% 

308 


158,150] 

[184,124] 

E4 

106.75 

3 

- 0.002°{7") 

- 0.002°(6.3") 

1.23 

-7.87 

-6.1, -7.0 

NGC4552 (M89) 

10.78 

15.7 

-60% 

456 


204,252] 

[213, 243] 

SO.l.(O) 

85.93 

3.6 

- 0.002®(7.4") 

- 0.002®(6.4") 

1.25 

-8.02 

-6.3. -7.2 


Column 


lagnitude from |Cd^ et al.| |2004j. 

^Luminosity distance of the galaxy (from NED), 
fraction of the D 25 area of each galaxy covered by the ACSVCS ob8ervations^^__^^^^ 
total number of GC candidates for ejich galaxy according to the catalog of|jordan et al. || 

® number of GCs in the red/blue classes 
^number of GCs in the high-L/low-L classes 

® morphological classification of the host galm^^omlFerrMese_et al. 

^ effective radius re of the host galaxy from|^rrMese^^^^"^^^™~” 

‘ D 25 expressed in units of the effective radius re 
* size of pixels in the residual maps along the R.A. axis 
size of the pixels in the residual maps along the Dec. axis 
” threshold g — z color used to define the red and blue GC classes 
° threshold Mg absolute magnitude used to separate high-L and low-L GC classes 

P Limiting absolute magnitude of the GC catalog in the g filter, based on the completeness estimated by|Jord^^^e^^a^|200^ using simulations. The limiting magnitudes are shown for two 
radii rh = 0.012S'' and =0.1283^', assuming an intermediate background level (f>g=0.4023 mag/arcsec^). 

The fraction of area covered by observations and the number of GCs in brackets have been obtained from the imaging data and GCs catalog presented byjstrade^^^ ^aj|201^, consisting i 
pointing plus five complementary HST pointings of the outskirts of the galaxy obtained through a joint C/iandro-HST proposal (P.I. Fabbiano). 

The fraction of area covered by the observations and the number of GCs in brackets have been obtained from the imaging data and GCs catalog used byjBlo^^e^ ^alj|2012^, consisting i 
pointing plus six complementary HST pointings of the outskirts of the galaxy (P.I. Sivakoff). 


values of the half-light 
n the central ACSVCS 
n the central ACSVCS 


luminosity classes produce somewhat redundant samples. 
For instance, the overlap between red and high-L classes 
varies between 56% to 87% according to the host galaxy. 
For this reason, similarities in the residual maps of these 
classes are to be expected. We present the results for 
both color and luminosity classes for the sake of com¬ 
pleteness. 


3. METHOD 


Following 

D’Abrusco et al. 

(2013) (see also D’Abrusco 

et al. 2014al 

□), tor each sampl 

e of GCs, we generated den- 


of Dressier (1980). These maps are based on the local dis- 
tribution of GCs, i.e. on the distance of the A-th closest 
GCs from the points where the density is estimated. The 
density is evaluated in each knot of a regular grid cov¬ 
ering the region of the sky where GCs are found, and is 
given by: 


Dk 


K 

Aoid-K) 


( 1 ) 


where K is the index of the nearest neighbor used to 
calculate the density, and A^idK) = n ■ is the area 
of the circle with radius equal to the distance dK of the 
A-th nearest neighbor. The uncertainty on the KNN 
density scales with the square root of K, so that the 
relative fractional error is inversely proportional to the 
square root of A. For this reason, the fractional accu¬ 
racy of the method increases with increasing A at the 
expense of spatial resolution. In this paper, we recon¬ 
struct the spatial distribution of GCs using values of A 
ranging from 2 to 10. The size of the spatial structures 
depends on a combination of the value of A and the den¬ 
sity of sources. Therefore, the optimal value for A de¬ 
pends on how extended and densely populated are these 
structures. In our experiments we consider a range of 
A values that are roughly proportional to the average 
density of sources of the samples for which the density 
maps are reconstructed. Moreover, very large values of 
A would smooth over spatially smaller structures thus 
losing spatial information. 

For all the GC distributions, we produce residual maps 


by subtracting smooth azimuthally homogeneous density 
maps from the observed density maps. These smooth 
maps were generated to have the radial number density 
profiles of the observed GCs. Using MonteCarlo simu¬ 
lations, we then estimate the statistical significance of 
each pixel in the residual map in terms of deviation from 
the mean value of the simulated density distribution in 
the same pixels (after making sure that such distribution 
of simulated densities can be modeled with a Gaussian). 
Finally, we estimate the statistical significance of spa¬ 
tial structures containing multiple pixels, by taking into 
account the spatial features of the structure, i.e. the 
number of member pixels and the peculiar shape of the 
pixel aggregation. 

The density maps for all galaxies in our sample of 
galaxies were evaluated on a regular spatial grid with 
spacing varying between A(R.A) ~ 0.002° (^ 6") and 
A(R.A) ~ 0.005° (~ 18") along the R.A. axis, and Z\(Dec) 
-0.002° (-6") and ZM;R.A) -0.005° (- 19") along the 
Dec. axis. The sizes of pixels for each target (columns 
(i) and (1) of Table 1) were chosen to produce a similar 
number of pixels in the grid covering the regions occu¬ 
pied by the GCs in all the galaxies investigated. This 
strategy permits to meaningfully compare the statisti¬ 
cal reliability and the area measured in pixels of each 
spatial feature for different galaxies without taking into 
account the different physical scales of the host galax¬ 
ies. For each galaxy, we have also compared the residual 
maps obtained by randomly varying the pixel sizes along 
both axes within a ±30% interval centered on the values 
shown in Table 1, to rule out systematic effects on the 
determination of the spatial features. The density in the 
boundary pixels was weighted according to the fraction 
of the pixel within the observed region. 

In Appendix A we compare the KNN method with 
other density reconstruction techniques, and show how 
it is optimal for our analysis. 

4. STRUCTURES IN THE SPATIAL DISTRIBUTION OF GCS 

For each galaxy, we have performed experiments for the 
entire GC sample and for the red/blue and low-L/high- 
L classes. Each experiment involved 50,000 simulations. 
Table 2 shows the percentage of simulated distributions 
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of GCs with “extreme” number of pixels, i.e. with val¬ 
ues exceeding those of the 90% of pixels in the observed 
residual map for K = {7, 8, 9}, for each galaxy. The K 
values for which these fractions are small (< 1%) corre¬ 
spond to the spatial scale at which the observed struc¬ 
tures become statistically significant. The parameter K 
is a measure of the expected spatial scale of the inves¬ 
tigated structures and of the density contrast of these 
structures over the average local density. Higher val¬ 
ues of K are more suitable to detect structures located 
within the £>25 of the galaxies, where the overall density 
of GCs is larger, since only high-contrast structures can 
be reliably detected over this high-density background. 
Smaller K values, on the other hand, are more apt at 
detecting structures in regions where the total number 
of GCs is smaller, as in the outskirts of the host galaxies. 
Since the data used in this paper all cover the central re¬ 
gions of the galaxies, we have focused on large values of 
K. For the six galaxies NGC4472, NGC4486, NGC4649, 
NGC4406, NGC4526 and NGC4621, we find significant 
inhomogeneities for K = [9, 9, 8, 9,8,9], respectively. For 
NGC4382 and NGC4552 we find marginally significant 
inhomogeneity for K = 9. The GC systems of NGC4374 
and NGC4365 are not consistent with an overall inhomo¬ 
geneous distribution for the range of K values explored 
(see the non-negligible fraction of simulated distributions 
with similar overall number of “extreme” pixels in col¬ 
umn (1) in Table 2), however, even in these cases lo¬ 
calized statistically significant structures can still be ob¬ 
served. 

Based on the residual maps of the entire GC system 
evaluated for the optimal K values of the six galaxies 
NGC4472, NGC4486, NGC4649, NGC4406, NGC4526 
and NGC4621 and ioi K = 9 for the remaining four galax¬ 
ies, we have considered as “over-density structures” all 
the groups of contiguous over-dense pixels composed of 
at least one pixel with significance larger than 4a. The 
residual maps are shown in Figures through We 
calculate the total significance of a GG structure by tak¬ 
ing into account the average significance of the enclosed 
pixels, its size and the geometry. The total significance 
of a GC structure expresses the probability of observing 
a similarly shaped GC structure with comparable size 
and average significance in a simulated GC density map. 
More details o n the estimation of the total significance 
can be found in|D’Abrusco et al.| (|2013|). We have classi¬ 
fied the structures based on their total significance, size 
in pixels and morphology. 

The three classes of over-density structures are defined 
as follows: small structures: groups of 4 or more ad¬ 
jacent pixels with total significance > 4a, or groups of 
one, two or three adjacent pixels with total significance 
> 5a; medium structures: groups of 5 to 19 adjacent 
pixels with with total significance >4a; large structures, 
defined as groups of 20 or more adjacent pixels whose 
total significance is > 4a. All the remaining groups of 
over-dense pixels not meeting the minimal requirements 
were discarded. We found 24 medium and 14 “large” 
GC over-density structures, while the total number of 
structures is 229. In the following analysis, we will focus 
on the medium and large over-density structures because 
the small structures could be produced by stochastic fluc¬ 
tuations in the number of GCs in just one pixel. Table 3 
displays the number of all the over-density structures 


found in the residual maps for the K values in Table 2, 
split according to their class, and the size (in number 
of pixels) for the medium and large structures. Table 3 
gives also the fractions of pixels of the residual maps con¬ 
tained in all the structures and only in the large ones. 
Table 4 contains the basic properties of the 38 medium 
and large over-density structures, including the area of 
the structures in pixels, the total and excess number of 
GCs located within each structure. In the following, we 
will describe the main density features observed in each 
galaxy, listed in order of decreasing brightness. 

4.1. NGC4472 

The K = 9 residual map of the spatial distribution 
of all GCs in NGC4472 (Figure upper-mid panel) 

is characterized by a significant 6.5a) “curly” over¬ 
density structure located N of the center of the galaxy 
(structure Al). Al lies almost entirely outside the II 25 
2.9 Tg) but it bends inward at its eastern end, radi¬ 
ally following the major axis to one fourth of its length; 
it is more prominent in the red GCs (Figure J h upper 
right panel) than in the blue GCs (Figure lower left 
panel). The high- and low-L residual maps (Figure 
lower mid and right panels) suggest spatial segregation, 
with low-L GCs mostly located in the radial section of 
the structure, and high luminosity GCs located outside 
of the 1125, not showing a clear radial dependence. We 
can rule out dynamical friction as cause of the observed 
spatial segregation since, in that case, we would have ob¬ 
served separation between the high-L, more massive GCs 
and the low-L, less massive GCs, with the former closer 
to center of the galaxy and the latter at larger galacto- 
centric distances. 

The smaller distinct over-density structure A4 (~4a) 
in the N-W corner could be an extension of Al. The 
second and third more significant density features after 
Al, labeled as A3 (~4.1a) and A2 (~4.3a) in the upper- 
mid panel of Figure are located N and W of the center 
of the galaxy. They are visible also in the residual maps 
for the two color classes and the high-L GCs. The first of 
these structures has a radial shape pointing N-W along 
the major axis. 


4.2. NGC4486 (M87) 

The highly inhomogeneous 2D distribution of all GCs 
in NGC4486 (Table 4, Figure]^, is dominated by a large 
over-density structure (Bl, ^ 5.8a), located in the N- 
E corner of the covered area (Figure upper left and 
mid panels). Bl is also visible, with sli^tly smaller sig- 
nihcance, in the color and luminosity GC classes (Fig¬ 
ure^ upper-right and lower panels). The second most 
signmcant over-density structure, in the S-E corner (B2, 
^ 4.8a), is also observed in red and low-L GCs. Mul¬ 
tiple over-density structures (B5, B6, B3, B7 and B4 
in Eigure upper mid panel) follow the major axis of 
NGC4486, turning to the N-W corner of the field, sug¬ 
gesting the presence of a large scale radial structure. 
Four structures (> 4.2a) are also visible in high-L GCs 
(Figure!^ lower mid panel). Three of these, in the N-E 
co rner of the field, correspon d to Bl, B3 and B6. 

|Romanowsky et al.| (2012) discussed two kinematically 
selected GG substructures m M87, one associated to the 
outer-halo stellar filament N-W of the galaxy, and the 















TABLE 2 

Fractions of simulated residual maps with number of “extreme” pixels exceeding the thresholds. 
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Galaxy 


All GCs 

(a) 

Red GCs 

(b) 

Blue GCs 
(c) 

High-L GCs 

(d) 

Low-L GCs 
(e) 

NGC4472 (M49) 

K = 7 

12.1%(3.1%) 

13.7%(T0%) 

9.2%(3.8%) 

9.9%(4.0%) 

12.4%(2.1%) 


K = 8 

4.2%(0.7%) 

8.6%(0.5%) 

6.6%(0.5%) 

4.3%(0.2%) 

3.9%(1%) 


K = 9 

1.1% (0%) 

2.3%(0%) 

0.9%(0%) 

0.1%(0%) 

0.3%(0.1%) 

NGC4486 (M87) 

K = 7 

10%(2.3%) 

12%(4.1%) 

8.8%(3.5%) 

11%(4.4%) 

7.5%(4.4%) 


K = 8 

3.4%(0.7%) 

8.7%(1.9%) 

6.6%(1.5%) 

9.1%(2.1%) 

1.5%(0.9%) 


K = 9 

0.6% (0%) 

1.4%(0.2%) 

2.9%(0.2%) 

4.3%(0.1%) 

0.2%(0%) 

NGC4649 (M60) 

K = 7 

3.3%(0.2%) 

0.2%(0.3%) 

1.2%(0%) 

3.6%(0.5%) 

1.9%(0.7%) 


K = 8 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 

0.1%(0%) 


K = 9 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 

NGC4406 (M86) 

K = 7 

4.1%(0.5%) 

5.2%(0.9%) 

9.3%(3.1%) 

3.1%(0.4%) 

1.9%(0.4%) 


K = 8 

1.4%(0%) 

1.7%(0%) 

5.5%(1.2%) 

0.2%(0%) 

0%(0%) 


K = 9 

0%(0%) 

0.1%(0%) 

1.4%(0.1%) 

0%(0%) 

0%(0%) 

NGC4382 (M85) 

K = 7 

13.4%(5.4%) 

12%(3.3%) 

18.9%(9.1%) 

7.9%(2.1%) 

9.3%(3.1%) 


K = 8 

8.9%(2%) 

8.6%(0.8%) 

11.1%(3.7%) 

4%(0.1%) 

4.5%(1.2%) 


K = 9 

3.3%(0.2%) 

2.2%(0%) 

9.2%(0.9%) 

0.2%(0%) 

0.6%(0.1%) 

NGC4374 (M84) 

K = 7 

24.4%(11.2%) 

11.6%(7.8%) 

14.5%(5.5%) 

19.9%(8.8%) 

13.5%(8.1%) 


K = 8 

21%(9.7%) 

9%(5.1%) 

3.9%(1.1%) 

12.3%(4.6%) 

9.7%(4.4%) 


K = 9 

12.5%(2.9%) 

4.9%(0.8%) 

2.1%(0.3%) 

7.8%(1.9%) 

5.4%(3.2%) 

NGC4526 

K = 7 

0.1% (0%) 

0%(0%) 

0%(0%) 

0%(0%) 

0.3%(0%) 


K = 8 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 


K = 9 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 

0%(0%) 

NGC4365 

K = 7 

36.6%(25.7%) 

30.8%(19.5%) 

21.3%(14.4%) 

18.3%(10.4%) 

40.6%(31.1%) 


K = 8 

31.4%(23.1%) 

28.2%(14.5%) 

15.6%(9.1%) 

12.7%(7.6%) 

26.9%(19.1%) 


K = 9 

24.3%(18.9%) 

17.6%(8.3%) 

5.9%(0.9%) 

5.7%(0.5%) 

22.2%(12%) 

NGC4621 (M59) 

K = 7 

8.4%(0.4%) 

3.9%(0.1%) 

4.8%(0.5%) 

6.5%(1.7%) 

1.9%(0%) 


K = 8 

1.1% (0%) 

0.5%(0%) 

2%(0%) 

4.3%(0.2%) 

0.3%(0%) 


K = 9 

0%(0%) 

0%(0%) 

0.1%(0%) 

0.2%(0%) 

0%(0%) 

NGC4552 (M89) 

K = 7 

12.9%(7.6%) 

13.6%(6.5%) 

7.7%(1.9%) 

6.4% (2%) 

13.7%(8.8%) 


K = 8 

6.5%(3.2%) 

9.5%(4.8%) 

3.6%(0.5%) 

3.2%(0.4%) 

9.5%(3.5.%) 


K = 9 

1.9%(0.2%) 

3.6%(1.3%) 

1.2%(0.1%) 

0.9% (0%) 

2.7%(0.8%) 


Column description: 

fraction of simulations obtained with K = {7, 8, 9} with number of extreme pixels with density exceeding 
the 90-th percentile of the observed pixel residual distribution) for the entire sample of GCs. 
as in column (a) for red GCs only, 
as in column (a) for blue GCs only, 
as in column (a) for high-L GCs only. 

® as in column (a) for low-L GCs only. 


other in the inner halo (the “shell”). This second struc¬ 
ture was interpreted as the signature of a disrupted in¬ 
falling system. Both phase-space structures are located 
outside of the region observed by A C SVC S. We have ap¬ 
plied our metho d to the sample of Romanowsky et al. 


(2012) (see also Strader et al. 
ity to detect these kinematica' 


2011 i to assess our abil- 

ly selected structures as 


over-densities in the GC spatial distribution (Figure 
lower panels). We sho uld keep in mind, h o wever , that 
the ACSVCS and the Romanowsky et al.| (|2012|) sam¬ 
ples have different level s of co mpleteness; in particular, 
the Romanowsky et al. (2012) sample is biased towards 
higher luminosity (mass) GCs. 

At very large radii, we find a high-significance overden¬ 
sity structure W of NGC4486, located where seven GC 
members of the shell (green points in lower, left panel 
in Figure]^ can be found. While the other members of 
the shell are distributed at the border of under-density 
areas (blue), they do not constitute a statistically signif¬ 


icant overdensity. A z oomed-in version of t he re sidual 
map obtained from the Romanowsky et alj(2012) GCs 
sample for K = 9 (lower mid panel, h'igure0 shows that 
a very significant (cr~ 10), extended GC structure is lo¬ 
cated N and N-W of NGC4486, mostly outside the D 2 B 
of the galaxy. A low-significance bridge of overdensity 
pixels connects this large GC structure to the complex 
of GCs over-densities observed in the ACSVCS residual 
map (Bl, B3, B4 and B7). The avera ge line-of-sight 
velocity of all Romanowsky et al. (2012) GCs enclosed 


within this large structure is consistent with the bulk ve- 
loci ty of the system within 2a. The average velocities of 
the Romanowsky et al. (2012) GCs located within the 
ACSVCS GC structures are consistent with the bulk ve- 
loc ity of the system, but the numbers of GCs are small 
(4 Romanowsky et al. 2012 GCs in Bl, 3 in B2, 2 in 
B3, 8 in B4, 4 in B5, 3 in B6 and 0 in B7). For this 
reason, we cannot rul e out a distinct kinematics, based 
on the Strader et al. (2011) sample only. A complete 
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TABLE 3 

Global properties of the over-density structures. 


Galaxy 

Ntot 

(a) 

^medium 

(b) 

large 

(c) 

% observed area in 
medium -|- large 
structures (d) 

% observed area in 
large structures 

(e) 

% GCs in 

all structures 

(f) 

% excess GCs in 
all structures 
(g) 

NGC4472 (M49) 

27 

3 

1 

UWo 

3 : 4 % 

IT?7% 

49 % 

NGC4486 (M87) 

46 

6 

1 

16.9% 

2.1% 

17.2% 

6 . 5 % 

NGC4649 (M60) 

28 

0 

1 

16.6% 

8.5% 

10.9% 

5 . 7 % 

NGC4406 (M86) 

12 

3 

1 

16.5% 

11.7% 

26.4% 

13.2% 

NGC4382 (M85) 

21 

2 

2 

14.6% 

7.2% 

22.3% 

9.9% 

NGC4374 (M84) 

22 

3 

1 

13.6% 

3.2% 

14.4% 

6.8% 

NGC4526 

4 

0 

2 

30.1% 

28.1% 

49.2% 

28.6% 

NGC4365 

38 

4 

1 

15.2% 

2.8% 

9.7% 

4.3% 

NGC4621 (M59) 

10 

0 

2 

18.4% 

16.2% 

35.4% 

16.4% 

NGC4552 (M89) 

21 

3 

2 

16.6% 

5.8% 

21.3% 

9.8% 


Column description: 

^ number of all GC over-density structures, 
number of medium GC over-density structures, 
number of large GC over-density structures 

fraction of observed area of the galaxy belonging to medium and large GC structures. 

® fraction of observed area of the galaxy belonging to the large GC over-density structures. 
^ percentage of GCs located within GC structures over the total number of GCs 
s percentage of excess GCs located within GC structures over the total number of GCs 



Right Ascension [deg] 


Right Ascension [deg] 




Right Ascension [deg] 



Right Ascension [deg] 


Fig. 1. — (Upper left) Scatterplot of the positions of the GCs detected in NGC4472 and used in this paper. Orange and blue points 
correspond to red and blue GCs, while solid and open symbols show high-luminosity and low-luminosity sources, respectively, defined using 
the magnitude threshold discussed in Section]^ (Upper middle) Residual map obtained for K = 9 of the whole sample of GCs in NGC4472. 
Blue and orange pixels show negative and positive excess, i.e. under-density and over-density residual. The intensity and symbols express 
the degree of significance of each pixel: light blue/orange, dark blue/orange, horizontal bars and crosses indicate over/under-density pixels 
with significance smaller or equal than 1, smaller or equal than 2, smaller or equal than 3 and larger than 3, respectively. (Upper right) 
K = 9 residual map for red GCs. (Lower left) K = 9 residual map for blue GCs. (Lower middle) K = 9 residual map for high-luminosity 
GCs. (Lower right) K = 9 residu al map for low-luminosity GCs. In all plots, the ellipse shows the D 25 (~2.9re) isophote of the galaxy 
from [de Vaucouleurs et al.| ( |199T] |. The solid black line shows the footprint of the HST observations used to extract the GCs. 
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TABLE 4 

Main properties of the “medium” and “large” GCs over-density structures. 


Galaxy 

Structure 

(a) 

Class 

(b) 

Significance 

(c) 

Size [pixels] 

(d) 

Total GCs 

(e) 

Excess GCs 

(f) 

Physical size [kpc] 
(g) 

NGC4472 (M49) 

A1 

All, red, blue, high-L*, low-L"^ 

~ 6.5(7 

22 

25 

8.1 

14.4 


A2 

All, red, blue, high-L, 

~4.3<7 

9 

46 

23.2 

5.3 


A3 

All, red, blue, high-L, low-L 

~4.1<7 

10 

47 

16.6 

4.6 


A4 

All, red, high-L 

~4<7 

7 

6 

2.1 

1.1 

NGC4486 (M87) 

B1 

All, red, blue, high-L, low-L 

~ 5.8^7 

18 

23 

11.3 

5.7 


B2 

All, red, blue, low-L 

~4.8<7 

20 

31 

12.1 

1.2 


B3 

All, red, blue, high-L 

~ 4.7^7 

9 

15 

6.8 

1.4 


B4 

All, red, blue, low-L 

~4.5<7 

9 

83 

38.7 

5.8 


B5 

All, red, low-L 

~4.5<7 

6 

28 

8.8 

1.6 


B6 

All, red, high-L, low-L 

~4.3<7 

7 

25 

12.8 

4.6 


B7 

All, red 

~4.1(7 

18 

18 

5.6 

1 

NGC4649 (M60) 

Cl 

All, red, high-L*, low-L* 

~10f7 

60 

160 

80.3 

29 

NGC4406 (M86) 

D1 

All, red, blue, high-L, low-L 

~ 15(7 

91 

76 

38.8 

14.3 


D2 

All, blue, high-L 

~4.5<7 

7 

4 

1.2 

2.4 


D3 

All, red, low-L 

~4.3<7 

8 

5 

2.4 

4.6 


D4 

All, red, high-L 

~4(7 

8 

4 

2 

5.1 

NGC4382 (M85) 

El 

All, red, blue, high-L, low-L 

~ 8.5^7 

61 

36 

18.6 

7.9 


E2 

All, red, blue, high-L*, low-L* 

~ 6.8^7 

22 

26 

11.8 

6.4 


E3 

All, red, high-L 

~ 6.3(7 

12 

20 

3.4 

3.1 


E4 

All, red, blue, high-L, low-L 

~ 5.7^7 

19 

19 

9.1 

6.2 

NGC4374 (M84) 

FI 

All, blue, high-L 

~ 6.5(7 

30 

18 

9.2 

6.4 


F2 

All, red, blue, high-L, low-L 

~ 6.3(7 

19 

9 

3.2 

3.5 


F3 

All, blue, high-L 

~ 5.5(7 

17 

23 

6.6 

3.6 


F4 

All, red, high-L, low-L 

~4.5<7 

11 

7 

1.4 

2.7 

NGC4526 

G1 

All, red, blue, high-L, low-L 

> 20(7 

185 

63 

37.9 

9.3 


G2 

All, red, blue, high-L, low-L 

> 20(7 

88 

53 

28.8 

16.7 

NGC4365 

HI 

All, red*, blue* 

~ 5.8^7 

22 

72 

29.5 

23.2 


H2 

All, red, blue, high-L, low-L 

~ 5.1(7 

10 

45 

5.6 

18.7 


H3 

All, blue, high-L, low-L 

~4.6(7 

5 

24 

6.2 

11.2 


H4 

All, blue, high-L, low-L 

~4.5<7 

6 

28 

18.4 

4.3 


H5 

All, red, blue, low-L 

~4.1(7 

6 

24 

12.5 

4.2 

NGC4621 (M59) 

11 

All, red, blue, high-L, low-L 

> 15(7 

71 

42 

23.4 

10.9 


12 

All, red, blue, high-L, low-L 

> 15(7 

85 

56 

21.3 

7.7 

NGC4552 (M89) 

LI 

All, blue, high-L 

~ 8.7^7 

30 

7 

4.1 

5 


L2 

All, red, blue, low-L 

~ 7.3^7 

22 

11 

6.4 

4.2 


L3 

All, red, blue, high-L, low-L 

~ 6.4(7 

18 

26 

10.8 

4.3 


L4 

All, red, blue, high-L, low-L 

~ 5.7^7 

10 

23 

11.2 

3.3 


L5 

All, red, high-L*, low-L* 

~4.8<7 

17 

16 

8.5 

4.2 


Column description: 

^ unique code associated to the GCs structure 

^ significance of each structure in the residual maps generated for all GCs, the color and luminosity GC classes: boldface text is used for the 
GC classes where the structure is clearly visible and has large statistical significance, plain text is used for GC classes where the structure is 
still detected with lower significance, while the GC classes where the structure is not detected - or detected in a distinctly different shape and 
size - are not reported. Asterisks indicate cases where the residual maps of either the color or luminosity GC classes suggest the presence of 
spatial segregation 

^ statistical significance of the structure 
^ size of the structures measured in pixels 
® total number of GCs within the structures 

^number of GCs within the structures in excess to the smooth distribution of GCs used to evaluate the residual maps 
s approximate maximum physical size of the GC structures [kpc] 



photometric coverage of M87, at least including the area 
within its D 25 and a deeper, more homogeneous spec¬ 
troscopic campaign that can match the ACSVCS sample 
are needed to fully map its GC streamers and conhrm 
their nature in the phase-space. 


NGC4486 shows dust filaments (Ferrarese et al. 20061 


and has been detected by|di Serego Alighieri et al.|(2013l 
in the far infrared (FIR) images at 250 /rm obtaiiied by 
the Herschel Virgo Cluster Survey. The FIR emission 
indicates the presence of diffuse dust that may affect 
the detection of GCs and, in turn, our reconstruction of 
their spatial distribution. The location of the observed 
dust structures is not associated with significant under- 
or over-density regions in our residual GC maps, sug¬ 
gesting an overall negligible effect of the dust on the the 
reconstruction of the GC spatial distribution obtained 
with our method. 


4.3. NGC4649 

The residual map (AT = 8) of t he spatial distribution of 
all GCs detected in NGC4649 (|D’Abrusco et al.||2014al 
(Figure ^ upper mid panel) is dominated by a highly 
signihcant (~i0cr), spatially extended (~50 pixels) over¬ 
density structure (Cl) in the E side of the galaxy. This 
feature, less prominent in red GCs, originates close to 
the major axis, following it out to the D 25 (~4.5 Vg) be¬ 
fore bending northward; high-L and low-L GCs appear 
spatially segregated in Cl (Figure!^ lower mid and right 
panels), with the former mostly Ideated along the ma¬ 
jor axis and the latter following the D 25 ellipse. The 
blue GCs residual map (Figure [r lower left panel) also 
shows few more marginally significant over-density struc- 
tu res in the S and W sides of the galaxies. As discussed 


D’Abrusco et al. (2014a), the GC over-density struc¬ 


ture is not associated to any feature in th e diffuse stel¬ 


lar lig ht. Recent kinematic measurements (Arnold et al. 
2014) of the diffuse stellar light have revealed disk-like 


outer rotation, which suggests a major dry merger un¬ 
dergone by a massive lenticular galaxy progenitor. No 
obvious correlation between the GCs structure and the 
radial velocit ies of a sample of GCs obs erved spectro- 
scopically by iLee et al.l (|2008D was found (iD’Abrusco et 
aL]|2()14a[ ). 


4 .4. NGC 44 O 6 

The K = 9 residual map of all GCs in NGC4406 (Fig¬ 
ure]^ upper mid panel) is dominated by a large, spatially 
coherent structure (Dl, ^ IScr) located N-E of the cen¬ 
ter of the galaxy. At small galactocentric distances, Dl 
follows the radial direction, but it widens as it reaches 
the boundary of the held. Dl is clearly visible in both 
color and luminosity GC classes (Figure upper right 
and lower panels). While Dl in red GCs is mostly lo¬ 
cated along the radial direction, in blue GCs it is more 
azimuthally extended. The second most signihcant resid¬ 
ual structure (D2, ^ 4.5a) overlaps the S-E end of the 
major axis. D3 and D4 (^ 4.3tT and ^ 4tT respectively) 
are visible in the whole sample of GCs, to the W of the 
major axis. D3 is more evident in the red (Figure]^ up¬ 
per right panel) and low-L (Figure ^ lower left panel) 
GCs, indicating the presence of a cluster of metal-rich, 
low-luminosity GCs at galactocentric radius ~1.4'. 

A sub-a rcsecond dusty disk has been observed in 
NGC4406 (Ferrarese et al. 2006 di Serego Alighieri et 


al.|2013D . Since the spatial scale of the structure that we 


detect in the GC distribution is signihcantly larger that 
few arcseconds, the presence of the dust feature should 
not affect our results. 


4.5. NGG4382 

The K = 9 spatial distribution of all GCs in NGC4382, 
(Figure^ upper mid panel), is characterized by the pres¬ 
ence of two large, signihcant over-density structures. El 
(> lOtr), located in the S-W corner of the HST held, is 
radial in the entire GC sample and is visible in all GC 
color classes (Figurej^ upper right and lower left panels). 
E2 (^6.8cr), S of the center of the galaxy, crosses the ma¬ 
jor axis at ^ 0.8' from the center. E2 and El could be 
distinct sections of a single elongated structure. A sec¬ 
ond group of signihcant and spatially coherent structures 
is located N-W; these are clearly visibile in the residual 
maps of all and red GCs: E3 (~6.3cr) and E4 (~5.7 (t). 
E4 looks similar to El, and extends radially from the 
center to the N-W corner of the held (Figure upper 
mid panel). Another property of the spatial distribution 
of all GCs in NGC4382 is the very signihcant (^ 6.5cr) 
E-W anisotropy, with ^ 68% of all GCs located E of the 
major axis of the galaxy. Since the ACSVCS image is 
not centered on NGC4382 center, it is not surprising tat 
the asymmetry between the E and W equal area halves 
of the ACSVCS held is slightly more signihcant, with 
^72% of the total GCs located in the western half of the 
held. 


4.6. NGG 4374 

This galaxy has only one large over-density structure, 
FI (~6.5 (t), located in the S-W corner of the ACSVCS 
data held (Figure]^ upper mid panel). FI is not visible in 
the red and low-luminosity GC classes (Figure]^ upper 
right and lower mid panels). The smaller over-density 
structures F2 (~6.3tT), F3 (^S.Scr) and F4 (^4.5cr) are 
roughly located along the major axis. F3 lies nearby FI, 
suggesting that they may be physically connected. 


NGC4374 is known for the 

presence of abundant 

dust ( 

di Serego Alighieri et al. 

2013 

). In particular. 

two parallel dust lanes, located c 

ose to the center ot the 


galaxy and slightly N of it along the E-W direction (Fer¬ 
rarese et al. 2006), extend for ~ 5" and ~ 14", respec- 
tively. While the presence of the dust lane may have 
affected the detection of GCs close to the position of F4, 
the GC structures FI, F2 and F3 are located too far from 
the dust to be signihcantly affected. 

4.7. NGG4526 

We detected (Figure upper mid panel) two strong 
spatially extended over-density structures located W 
(Cl) and E (G2) of the center of the galaxy, both > lOcr 
signihcant and clearly visible in all the classes (Figure]^ 
upper right and lower panels). Cl, at small galactocen¬ 
tric distances (< 1'), roughly follows the W section of the 
major axis of the galaxy, and it extends towards N as it 
approaches the boundary of the held covered by the HST 
data, at galactocentric distance between ~1' and ~1.6'. 
G2 overlaps the E section of the major axis. 


Since NGC4526 is a lenticular (S0_3_(6)) galaxy (Fer- 


rarese et al. 

2006 

), the two structures Cl and C2 could 

be enhancec 

by projection ettects, assuming a smooth 
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Right Ascension [deg] 


Blue GCs 
K=fi 



Right Ascension [deg] 



Right Ascension [deg] 




Right Ascension [deg] 





Fig. 2.— Upper and mid panels: scatterplot of the position of the GCs in NGC4486 and residual maps obtained for K = 9 for the whole 
sample of GCs and the two color and luminosity GC classes. The D 25 elliptical isophote (~2.9re) is not shown because it falls outside of 
the plot borders. The major axis of the galaxy is shown by the vertical red line. Refe r to Figure for a description of the other p anels. 
Lower panels: residual map for K = 9 for the sample of GCs observed spectroscopically ([Strader et aL|2011| [Romanowsky et al.|2012|l. The 
GCs belonging to the shell selected kinematically are shown as green diamonds. On the rigut, a zoomed-in version ot tne same map with 
the region covered by the ACSVCS data and the GC structures discussed in Section |4.2| highlighted. 


three-dimensional distribution of GCs that follows the 
galaxy light. In order to check whether projection ef¬ 
fect may bias our results, we have produced additional 
residual maps using a modified model for the simulated 
spatial distribution of GCs. We add to the smooth back¬ 
ground component with isotropic azimuthal distribution 
(the halo) an azimuthally anisotropic component aligned 
along the major axis, spanning 15 deg above and below 
the major axis measured at the ZI 25 ) to mimic the pres¬ 
ence of a population of GCs associated to the disk. We 
have performed different experiments where the fractions 
of simulated GCs associated to either component varies 
between 30% and 70%. Both G1 and G2 are still visi¬ 


ble, although with decreasing significance as the fraction 
of simulated GCs in the disk component grows. In the 
worst case scenario, the significance of G1 is ~ 7cr and 
the significance of G2 is ^ 3.5(T.Based on the results, we 
can rule out that both G1 and G2 structures are entirely 
spurious. 


Moreover, NGC4526 contains abundant du st (Fer- 
rarese et al.|2006 di Serego Alighieri et al. ]2M^ , mostly 


located in a large-scale (15.6") disk structure aligned 
along the major axis of the galaxy. The non-negligible 
dust obscuration may have affected the detection effi¬ 
ciency of GCs and, in turn, the characterization of their 
spatial structure. Spectroscopic measurements could 
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Fig. 3.— Scatterplot of the position of the GCs in NGC4649 and residual maps obtained for K = 8 for the whole sample o f GCs and 
the twq_color and luminosity classes. These plots are slightly modified versions of the plots shown by[D’Abrusco et al.[ ([2014a]. Refer to 
Figure]^ for a description of each panel. The smaller ellipse represents the area o ccupied by NGC46y . a less massive spiral galaxy that is 
thought to be interacting with NGG4649. This area was not observed with ACS ([Strader et al.|2012[|. 


confirm the nature of G1 and G2 as physically distinct 
groups of GCs by detecting them also in the phase-space. 

4.8. NGCi365 

The spatial distribution of all GCs detected in 
NGC4365 (Figure upper mid panel) is characterized 
by the absence of Targe structures and by a marginally 
significant (2.8 (t) asymmetry in the distribution of GCs, 
which are more numerous to the E of the major axis of 
the galaxy. The most significant feature HI (~ 5.8cr), 
is close to the intersection of the £>25 S.&’e) ellipse 

with the major axis. HI extends towards the N-E corner 
of the field, where H4 4.5cr) can be observed. HI is 
visible in the residual maps for both red and blue GC 
classes (Eigurej^ upper right and lower left panels) with 
a strong indication of spatial segregation: red GCs in the 
W of the L structure and blue GCs in the E. H2 (~ 5.1 (t), 
slightly E of the intersection between the £>25 and the S 
end of the major axis, is more prominent in the high- 
luminosity and blue GCs (Figure!^ lower right and left 
panels). H3 (~4.6 (t) is located E m H2. 

Both H3 and H2 lie along the S-W direc tion, where 
the comp anion galaxy NGC4341 is located. [Bogdan ^ 
(20121 reported that the centers of NGG4365 and 
iNG(J4342 are connected by a stellar-light stream. This 
stream is an indication of tidal interaction between 
NGC4365 and NGC4342 and corresgonds to a stream 
of GCs detected by Blom et al. (20141 with kinematics 
consistent wi th the 1NGG4342 G C system (see Eigu res 10 
and 11 from Blom et al. 20141. Using data from Blom 


et al. (2014), we find a marginally significant (^ 2.2a) 
GC overdensity consistent with the GC. H2 and H3 are 
roughly located at the E end of the GC stream. How¬ 
ever, only one GC of the sample of 31 selected to have 
line-of-sight velocities compatible wi th the motion of t he 
NGC4342 GC system (Eigure 10 in Blom et al.||2014) is 
located within the region occupied by H2. The average 
line-of-sight vel ocity of the spect roscopic GCs from the 
list presented by Pota et al. (2013|) located within H2 and 
H3 is compatible with the systemic velocity of the GC 
system of NGC4365. H5 (^4.1cr), N of the center of the 
galaxy, across the £>25 ellipse is also visible (with vary¬ 
ing significance) in the other GC color and luminosity 
classes. 


4.9. NGC4621 

Similar to NGC4526, the spatial distribution of GCs in 
NGC4261 (Figure upper mid panel) shows two large 
and strong structures along the major axis. II (> 12 (t), 
in the S corner of the field observed, has similar signifi¬ 
cance and shape in the color and low-L classes (Figure!^ 
upper right and lower panels), but it appears positionally 
shifted in the high-L GCs (Figure lower mid panel), 
12 is similarly significant in all the oasses, but its shape 
and size change considerably. For example, in the red 
(Figure!^ upper right panel) and high-L GCs (Figure!^ 
lower mid panel), the section of 12 that extends from the 
major axis towards N-E is barely visible, while in the blue 
(Figurej^ lower left panel) and the low-L GCs (Figurej^ 
lower ri^t panel), it splits in three distinct components 
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Fig. 4.— Scatterplot of the position of the GCs in NGC4406 and residual maps obtained for K = 9 for the whole sample of GCs and the 
two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 


that occupy most of the area of the HST field N of the 
center of the galaxy. 

NGC4621 has large ellipticity (E4 iFerrarese et al. 

TT 1 n.11 'tT I I H ^ J m 


20061, and displays a small stellar disk extending ~ 10'' 


in the W half of the galaxy at galactocentric distance 
^ 1', suggests the existence of a long, spatially coherent 
elongated sequence of structures extending from the N-W 
corner of the field to the center, possibly encompassing 
also L2, L3 and L4. 


is neglig 

ible ( 

Ferrarese et al.||2006 

di Serego Alighieri et 

Ferrarese et al. 

(20061 report the presence of two thin 

al. 

2013 

I. In order to estimate the influence of the pro- dust filaments ext 

ending ^ 6" from the center towards 


jection of an elongated three-dimensional distribution of 
GCs on the detection of the GC structures II and 12 , we 
have used the same approach described in Section 14.7[ 
The worst case scenario significance of the structure 11 is 
^ 5.8(7, while 12 disappears as a single structure when the 
fraction of simulated GCs belonging to the component 
aligned with the major axis or the galaxy becomes larger 
than 55%. Based on these results, we can assume that II 
is not entirely produced by projection effects while the 
real nature of 12 can only be determined by a more de¬ 
tailed modeling of the GC system, which is not possible 
with the data currently available. 

4.10. NGC4552 


NGC4552 (upper-mid panel of Figure 10) displays two 
large structures in the W side of the fielm LI (~8.7 (t), 
located in the S-W corner, extends parallel to the D 25 
S.Gce) ellipse for ~ 0.9' and is prominent in blue and high- 
L GCs. L2 (^7.3(7), slightly S of the N-W corner of the 
field, within D 25 , is clearly visible in all classes, except 
for the high-L GCs (Figure [lOl lower mid panel). Two 
less significant structures L3~(^ 6.4cr) and L4 (^ 5.7cr) 
(visible in all classes of GCs), indicate the presence of a 
significant over-density close to the center of the galaxy. 
The fifth most significant residual structure L5 (~4.8cr) 


N-NW direction and of patches of dust within a radial 
distance of 11". While these dusty regions could influ¬ 
ence the reconstruction of the spatial distribution of GCs 
around the positions of these structures, their size is too 
small to affect them significantly. 

5. DISCUSSION 

The results presented in this paper confirm that lo¬ 
calized spatial structures in the 2D distributions of GC 
systems are common. In the ten brightest Virgo cluster 
elliptical galaxies we have found GC over-density struc¬ 
tures with a range of amplitudes, sizes, and shapes. We 
measure linear dimensions ranging from ~ 1 to ~ 20 kpc 
and azimuthal sizes ranging from ^ 10° to ^90°. We ob¬ 
serve GC structures with both simple, roughly circular 
shapes, and more complex forms. Below we discuss the 
characteristics of these structures, within the constraints 
provided by the existing limited observational coverage. 


Bonfini et aL|2012 

D’Abrusco et al.||2013 

2014a|b|) 

use these results 

to set simple observation-based 


straints to the size of the accreted satellite galaxies, that 
may have originated these structures. 

5.1. Location and shapes of GC over-density structures 
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Fig. 5. — Scatterplot of the position of the GCs in NGC4382 and residual maps obtained for K = 9 for the whole sample of GCs and the 
two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 


Most of the over-density structures are located between 
^ 1' and 2' (0.2-0.8 Ve) from the center of the host galaxy 
where all the galaxies in our sampl e are re asonably well 
covered by the ACSVCS (Figures [^. While the 
decline in the number of structures detected at smaller 
galactocentric radii may be real, we cannot exclude the 
effects of either incompleteness of the GC catalog in the 
central area of the galaxies, or the decreased ability of 
our method to select structures in compact regions with 
very high GC density. The decrease in the frequency of 
structures for radii larger than 2' is definitely caused by 
the decrease of the galactic area covered by the ACSVCS 
observations. For r>2.5' only NGC4472, NGC4649 and 
NGC4365 have adequate HST coverage. Figure 1^ illus¬ 
trates the large scatter in the fraction of areas covered 
for the different galaxies in our sample (at IrgS the frac¬ 
tional coverage varies between 0% and 100% per galaxy, 
with ^ 35% of the total area covered) and the lack of ho¬ 
mogeneous coverage for galactocentric radii larger than 
0.5r„. 


In Figure we show the spatial distribution of the 
over-density structures summarized in Table 3. The 
significance-weighted central coordinates of the struc¬ 
tures are plotted in the plane generated by their galacto¬ 
centric radius, expressed as projected angular distances 
and in units of the effective radii of each galaxy, and 
azimuthal angle, measured counter-clockwise from the 
N direction of the major axis of the galaxy. We also 
show the areas of the galactocentric radius vs azimuthal 
angle plane that can be accessed with the GCs sam¬ 
ples used in this paper as shaded colored background. 


From Figure |11[ it is evident that within the annu¬ 
lus ^ 1' - 2' wFere coverage is good for most galaxies, 
the smaller structures appear randomly azimuthally dis¬ 
tributed in the [0°,360°] interval, suggesting that some 
of them might be the result of stochastic fluctuations, 
as discussed in Section]^ The medium and large struc¬ 
tures (the largest size points in Figure [TI|) cluster along 
the major and minor axes < 30° mia 6ff > 330° or 
150° <0Ar < 210°). We explore this effect further in Fig¬ 
ure |13[ where we plot the medium and large structures 
(see section 13 and Table 4) in the galactocentric radius 
vs azimuthafangle plane. 

The most extended large structure along the radial di¬ 
rection is Cl, in NGC4649, which spans the interval of 
galactocentric distances [0.5, 2.8] Ve- Other large struc¬ 
tures that extend mostly along the radial direction are II 
and 12, in NGC4621, from 0.2 to 1.1 Vg and from 0.1 
Te to 1.2 Te respectively, and Cl and G2 in NGC4526. 
All these structures are located along the major axis of 
their host galaxy in the whole range of radial distances 
covered by the ACSVCS observations. The distribution 
in angular position of the “large” structures (Fig¬ 


ure 13) shows that 9 out of 14 of them are located near 
the major axis of the host galaxies. When the areas of 
the structures are considered, the larger structures are 
preferentially located along the major axes of the host 
galaxies, accounting for ~ 86% of the total “structure 
area”. This results appea rs consistent with the conclu¬ 
sion of [Wang et al.j (|2013|, reporting preferential align¬ 
ment of CC system with the major axis in Virgo galaxies. 
However, 5 large structures (Bl, Dl, FI, LI, L2) follow 
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Fig. 6.— Scatterplot of the position of the GCs in NGC4374 and residual maps obtained for K = 9 for the whole sample of GCs and the 
two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 

confirmed by Janowiecki et al. (2010|), which have ob¬ 
served a complex, extended shell system in the diffuse 
stellar light, suggesting radial accretion of a companion 
galaxy. 


the minor axis. This minor axis effect was not previ¬ 
ously reported. It will need to be considered when trying 
to explain these results with theoretical simulations.The 
morphology of the host galaxy may play a role in these 
preferential orientations: the GC structures Gl, G2, II 
and 12 observed in the galaxies NGG4526 and NGC4261 


(classified a s S0_3_(6) and E4, respectively, by Ferrarese 
et al.|2006 1 are almost perfectly aligned along the major 
axis ot their galaxies. 

The large coherent structures may be suggestive of the 
disruption and merging of a sizable companion galaxy. 
For at least NGG4649, stellar kinematics suggest a ma¬ 
jor dry (dissip ationless) merger, with a lenticular galaxy 
as progenitor ( Arnold et al.|2014 ). However, we also note 
that the large structures may be slightly more prominent 
for the red GC subsystems. These large structures are on 
average 1.6 times more significant in the red GCs resid¬ 
ual maps than the medium structures (~ 3.5(T/pixel for 
large structures to ~2.2(T/pixel for medium structures), 
suggestin g some dissipative merging-related GC forma¬ 
tion (e.g., Bekki et al.||2002 ). 


NGC4486 has regular isophotes but the kinematics of 
its GC syste m indicates that it is still accreting massive 
companions (Romanowsky et al.||2012); further evidence 
of stripping of low-luminosity dwarfs along radial orbits 
come from the structures in its diffuse stellar l ight at 
large galactocentric radii (Janowiecki et al.||2010). 

NGC4649 does not show significant structures in its 
stellar light distribution which might be associated to th e 


GC over-density structure Cl (D’Abrusco et al. 2014a|). 
However, kinematic evidence based on the Hlanetary 
Nebulae (PN) and GC systems, suggests that a. merger 
has contributed t o the assembly of this galaxy ( Teodor- 
escu et al. 20 11| Das et a l. 201 1 Coccato et aE 2013). 
In particular, Das et al. 1 201 1[) showed that PINs and 


GCs might belong t o separate dynamica l systems at radii 


5.2. Observational evidence of the assembly history 

For all the galaxies in our sample with the exception 
of NGC4374, NGC4526 and NGC4621, evidence of past 
merging or accretion of companions has been reported in 
the literature. 

NGC4472 shows a fairly regul ar diffuse stellar distri- 
bution with elliptical isophotes (Ferrarese et al. 20061 
within few re, but it has been long known to be i nteract- 


larger than 12 kpc. Arnold et al. (2014) used kinematics 
of the diffuse stellar light to discover disk-like outer rota¬ 
tion that might have originated from a massive lenticular 
galaxy progenitor. 

NGC4406 also has small-scale streamers that can be 


explained b y stripping of dwarf companions (Janowiecki 


et al. 2010); furthermore, observations of the cold dust 


emission m the surrounding area has confirmed that 
NGC4406 has stripped materials from th e nearby spi¬ 
ral ga laxy NGC4438 by tidal interaction (Gomez et al. 


ing with its comp anion dwarf galaxy UGC7636 (McNa¬ 
mara et al. 1994). Its active accretion history has been 


2010 ). 


T7GC4382 has boxy isophotes (Ferrarese et al. 2006) 
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Fig. 7.— Scatterplot of the position of the GCs in NGC4526 and residual maps obtained for K = 9 for the whole sample of GCs and the 
two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 


and presents abundant fine structure in its diffuse stel¬ 


lar light (Schweizer et al. 19901, suggesting that it has 
undergone a major merger in the past 4-7 Gyr. 
NGC4374, whose opt ical isophotes are very regular in 


the ACSVGS images (Ferra rese et al. 2006 1, shows a 
rapidly rotating central disk (Bower et a!i.|l998 ) but does 


not present evidence of kinematically decoup led compo¬ 


nents and i t l acks significant tidal structures (Janowiecki 


et al. 2010). Arnold et al. (2014) found that the outer 
region ot the galaxy show little to none rotation, in agree¬ 
ment with the expectations for a low rotator. 

NGC4526 is a lenticular galaxy hosting a large dusty 
disk that extends along the major axis. Its isophotes 
become clearl y boxy at ^ 10" befo re becoming disky at 
70" radius ( Ferrarese et al.|[2006 ). 


NGC4365 is known tor its kinematically decoupled core 
rotating along the major axis, which coul d be result of 
a ma jor merger involving a spiral galaxy ([Davies et al. 


2001). This galaxy is in the p rocess of strippin g GCs 


from the companion NGC4342 ( ]Blom et al.poid ). 


NGC4621 shows a thin faint stellar disk in the 
ACSVGS images, but further evidence of a recent ma¬ 
jor merger or accretion event is missing. 



accretion history involving either multiple small accre¬ 
tions or a major merger. 


5.3. GC structures in color and luminosity classes 


Color and luminosity differences in the GC over- 
dens^ structures are noticeable in several galaxies (Sec¬ 
tion S Table 4). While some of these differences may be 
due to statistical noise, some are significant and may pose 
constraints on the origin and evolution of the structures. 
Among the medium and large GC over-density struc¬ 
tures, spatial segregation has been detected for red/blue 
GC classes in HI (NGC4365), and for high-L/low-L GCs 
in Al (NGC4472), Cl (NGC4649), E2 (NGC4382) and 
L5 (NGC4552). We selected the GC structures with spa¬ 
tial segregation by comparing the significances of GC 
structures in each of the color (red/blue) and luminos¬ 
ity (high-L/low-L) classes. The structures for which the 
difference of significance was larger than Scr were visu¬ 
ally inspected to ascertain the presence of spatial segre¬ 
gation. In NGC4472, the luminosity segregation in the 
large structures Al suggests an incoming tangential GC 
stream that then turns into the radial direction. In this 
stream, the low-luminosity GCs are prevalent in the in¬ 
ner radial direction, while the more luminous GCs are 
mostly found in the section of the structure tangential 


to the D-jf,. This evidence suggests that dynamical fric- 

tion (Fall & Zhang||2001 

Binney & Tremaine 

2008) can 

be ruled out as the only mechanism responsi 

ale tor the 


observed segregation, as it would act more vigorously on 
the luminous, more massive clusters that should be ob¬ 
served at smaller galactocentric distances than the less 
massive, low-L GCs. Kinematic observations are needed 
to shed light on the origin of the spatially segregated 
high- and low-L components of Al. We stress that, lack¬ 
ing observational confirmation or a detailed explanation 



































































































































































15 



Right Ascension [deg] 



Right Ascension [deg] 




Right Ascension [deg] 


Right Ascension [deg] 



Right Ascension [deg] 



Right Ascension [deg] 


Fig. 8. — Scatterplot of the position of the GCs in NGC4365 and residual maps obtained for K = 9 for the whole sample of GCs and the 
two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 


of the physical mechanisms responsible for the observed 
properties of this GC structure, the two components in 
A1 could have originated in distinct events. 

The spatial segregation of high-L/low-L GCs observed 
in Cl (NGC4649) could be tentatively explained as a re¬ 
sult of differential dynamical evolution of the GC color 
classes formed during the maj or merger experienced by 
the host galaxy, as discussed by D’Abrusco et al. (2014a). 
Color differences in the distribution of the C(J over- 
density structures may be related to the prevalent GC 
population of merging satellite galaxies, if these are the 
results of dry mergers. Gas-rich mergers may increase the 
localized red GC population. We note that the highly in¬ 
complete spatial coverage of our sample (see Figure 12) 
does not allow us to assemble a full picture of either tEe 
shapes or the characteristics of these GC features. 


5.4. GC structures and local galaxy density 

Figure shows the location of the ten galaxies stud¬ 
ied in this paper in the region of the sky covered by the 
Virgo cluster. We compare their position with the local 
galaxy density, which we have evaluated from the cata- 
log of confirm ed and candidate Virgo galaxies of lBinggeH] 
et al. (1985). This catalog is complete to Bt ~ 18 mag 


and covers ~ 140 deg . The large GC structures detected 
in each galaxy, color-coded based on their membership, 
are shown in the insets of Figure [T^ for reference. We 
calculated the average fraction ofAhe observed galaxy 
areas belonging to medium and large GC structures for 
both higher- and lower-density regions of the Virgo clus¬ 
ter using the data shown in Table Only 2.1% of 


the total observed area of the three galaxies (NGC4406, 
NGC4486 and NGC4374) located within the inner green 
line in Figure |14[ which is the isodensity contour corre¬ 
sponding to 50% of the peak density of the Virgo cluster 
galaxies, is occupied by large GC structures (16.9% if 
both medium and large structures are considered). The 
fraction of the area of the remaining seven host galax¬ 
ies located outside of the 50% isodensity contour which 
is occupied by large structures grows to 8.5%, while it 
remains almost constant (16.6%) if both medium and 
large GC structures are considered. While these num¬ 
bers seem to indicate that galaxies in the highest-density 
regions of Virgo tend to have smaller GC structures than 
galaxies located elsewhere, we note that the uneven spa¬ 
tial coverage of ACSVCS is a source of uncertainty in 
these conclusions. For instance, if we consider only the 
five galaxies outside of the outer green line (the isoden¬ 
sity contour corresponding to the 33% of the peak den¬ 
sity), the fractions of area occupied by large and large 
plus medium GC structures become 2.8% and 15.2%, re¬ 
spectively. This inversion of the trend, though, is likely 
caused by the incomplete spatial coverage of the sample: 
NGC4649 and NGC4365, the only two galaxies almost 
completely observed within the II 25 , are both located 
outside of the 33% contour, and their presence “dilutes” 
the average fractional area occupied by large structures. 
Another possible source of bias is the uneven distribution 
of galaxy morphologies in higher- and lower-density re¬ 
gions: lenticulars and ellipticals with significant elliptic- 
ity are mostly located in the lower-density regions, while 
two of the three galaxies within the 50% isodensity con- 
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Fig. 9.— Scatterplot of the position of the GCs in NGC4621 and residual maps obtained for K = 9 for the whole sample of GCs and the 
two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 


tour (NGC4486 and NGC4374) show negligible to small 
ellipticity (EO and El, respectively). 

As discussed in Section |5.H the large GC structures 
are on average more significant 1.6 times) in the red 
GGs than in the general GC population. When split ac¬ 
cording to the environment density and normalized by 
the area of the structures, the large GC structures are 
even more significant in the residual maps generated by 
red GCs in galaxies in low-density regions than the large 
structures in medium and high-density areas of the clus¬ 
ter. In particular, the large GCs structures associated 
to host galaxies outside of the 33% of the peak density 
region are, on average, 5.3 times more significant than 
GC structures in all GCs, while the same excess proba¬ 
bility decreases to 3.5 and 1.2 for large structures in the 
medium and high-density regions. 

Our findings are in agreement with two different re¬ 
cent accretion histories for the galaxies in the higher- 
and lower-density regions of the cluster. As suggested by 
hydrodynamical sim ulations (Millennium simulation, see 


Springel et al. 20051 and observations of the l ocal clus¬ 


ter s of g alaxies (tor instance, for Coma cluster Smith et 


al. 20121, major merging events occurred earlier in the 


high-density regions at the center of the clusters. For 
instance, according to the Millennium simulation, an av¬ 
erage galaxy currently observed at small clustercentric 
distance joined a halo whose size is typical of the Virgo 
cluster 8.9 Gyr ago, corresponding to redshift ^ = 1.3. At 
a distance of 2 Mpc, the same galaxy reach ed a cluster¬ 


like e nvironment only 4.8 Gyr ago ( 2 ; = 0.5) (Smith et al. 


20121. Observationally, mergers in the distant past have 


depleted of similar-sized companion the environment of 
the largest galaxies, which are currently observed in the 
deepest potentials. Therefore, only recent accretions of 
smaller-than-average galaxies (satellites) can be expected 
in the densest areas of clusters. Instead, major mergers 
(or wet mergers with larger mass ratios) may have oc¬ 
curred in the more recent past in lower density areas, 
producing the large scale, spatially extended GC struc¬ 
tures mostly located along the major axis of the host 
galaxies, that we observe. 

5.5. Constraints on potential progenitors of the GC 
structures 


Th e ACDM mod el _ of hierarchical galaxy fo rma- 

tion (White & Rees 1978 Di Matteo et al. 2005) pre¬ 
dicts continuous galaxy evolution via merging and ac¬ 
cretion of satellites. Supporting evidence are provided 
by recent observations of the local Univer se, including 
the Sagittarius Stream (Ibata et al. 1994 ); Milky Way 


compani ons undergoing tidal disruption (Belokurov et 
|al.| | 2006|, and streams and dwa rf galaxies in the halo of 
M31 (jiVlcConnachie et al. 2009). The kinematics of the 
GC system of 1NG G4486 also suggest y five accretion of 
satellite galaxies (Romanowsky et al. 2012), as do the 
presence of GC streams in a f ew ellipticals (e. g.,|Stradeii 
^alj20lT Blom et al.|2012a). Other results (Rontini et 


al.||2012 ) and our pr evious work on the 2D GG distri bu¬ 


tions of E galaxies ([D’Abrusco et al.||2013[ |2014a|b|) as 
well as this work suggest that GG systems niay bear the 
imprint of these effects, as fos sil remnants of accretio n 
events of small dwarf galaxies ( Penarrubia et aL]|2009 ). 

Based on our results, we assume that the structures we 















































































































































17 



Right Ascension [deg] 


Right Ascension [deg] 



Right Ascension [deg] 



Right Ascension [deg] 



Fig. 10. — Scatterplot of the position of the GCs in NGC4552 and residual maps obtained for K = 9 for the whole sample of GCs and 
the two color and luminosity GG classes. Refer to Figure^ for a description of each panel. 


have detected in the 2D GC distributions of the ten most 
luminous Virgo ellipticals are the remnants of accreted 
galaxies, detected over the smooth, relaxed GC distri¬ 
bution of the main galaxy. A discussion of the caveats 
to this assumption can be found in Section |5.6| Using 
the excess number of GCs in the detected over-density 
structures, we can set lower limits on the masses of their 


parent galaxies, by using the relations of Harris et al.j 
(2013) between the number of GCs and the dynamical 
mass of elliptical galaxies. We stress that these esti¬ 
mates are lower limits, because GC s may be destroyed 


or reduced in these encounters (see Fall & Zhang 2001 
for an extensive review on the dyhamical evolution o: 
GCs). Moreover, detailed numerical simulations of the 
evolution of GCs in these scenarios are lacking. 

We considered all the GC medium and large over¬ 
density structures of the ten Virgo cluster galaxies in 
our sample (Table 5). Among these structures, we will 
focus, in particular, on the over-density structures that 
are clearly identifiable also in the residual maps obtained 
from the blue GCs, because metal-poor GCs could have 
been either produced early in the assembly history of the 
host galaxy or subsequently added thr ough accretion of 


low-mass, metal-poor sa tellite galaxies (Brodie & Strader 

These structures are A1 


2006 


Peng et al. 


2008|). These structures are Al, A3, 
B3, B4, HB, D2, E1,E2, E4, El, F2, F3, Gl, G2, 
HI, H2, H3, H4, H5, LI, L2, L3, L4 (the other structures 
for which the satellite properties have been calculated for 
comparison only are marked with an asterisk in Table 5). 

To estimate the number of GCs in each over-density 
structure, we used both the observed radial density pro¬ 


file Vexc,GC (which is a measure of the azimuthally av¬ 
eraged radial GC distribution in a ^laxy) and the min¬ 
imal radial density profile obtained by using 

the minimum value for each given azimuth (Table 5). 
These two different definitions of the density radial pro¬ 
files provide reasonable estimates of the minimum and 
maximum size of the GC systems in the structures, and 
were both used to evaluate the range of possible values 
of the basic properties of the accreted galaxies. We then 
corrected these estimates for incompleteness using the 
GC luminosity functions determined with the ACSVCS 
survey data for Virgo cluster galaxies by [Jordan et al.| 
(2007). The average additive corrections are ~ 0.6 and 
~U.9 GCs per density structure for the two radial density 
profiles used, respectively. We then applied the relations 


and VGCs = (refTe)i-29±0 03 ([Harris et 


al.|2013|) , where M^vn is the dynamical mass ([Wolf et al. 
^Vefinedas: - 


Mdyn = 


G 


( 2 ) 


which approximates the total baryonic mass of the galac¬ 
tic bulge within the effective radius where t he con- 
tribution of the dark-matter halo is negligible (Graves 


The minimal radial density (MRP) profile is defined as the 

density prohle obtained by setting the density in each radial bin 




cal terms^this assumption translates to assuming that all positive 


residuals in^h^^esiduafm^ ol (FJs cai^be_inter£rete^ajS_^ccreted 

Gu systems trom distinct satellite galaxies. ^ 
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Fig. 11.— Distribution of the over-density structures shown in Table 3 in the galactocentric radius vs azimuthal angle plane (measured 
counter-clockwise from the N direction of the major axis of the galaxy). In the left panel, the galactocentric radius i s expressed as projected 
angular distance, while in the right panel we use units of the effective radii re of each galaxy as estimated by |Ferrarese et al.| j2006|l 
(Table 1). In both panels, the shaded background regions represent the areas of the plots that can be accessed witn tne LiU sample used 
in this paper for each galaxy. These regions follow the same color-coding used for the points. In bot h pl ots, the large points represent the 
significance-weighted central positions of the medium and large GCs structures discussed in Section |5.1| and the gray and black marginal 
histograms show respectively the normalized distributions of the galactocentric radii of all GCs structures (gray) and large only GCs 
structures only (black). 



Fig. 12. — Fraction of the area of each galaxy as a function of 
the galactocentric distance in units of the effective radii re covered 
by the footprints of the HST ACS observations used to extract the 
catalogs of GCs employed in this paper. The black line represents 
the fraction of total area covered by the observations for the entire 
sample of galaxies. 


& Faber 2010); Ue is the central velocity dispersion of 
the galaxy derived from spectroscopic observations of the 


bulge available in the literature (see 

Harris et al. 

2013 

and references therein). Harris et al. 

(2013) found 

that 


TeCTf. strongly correlates with the size of the GC popula¬ 
tion. 


The results (Table 5) show that under our assump¬ 
tions, the average fraction of mass contributed by each 
accretion event (among the GC structures with no aster¬ 
isk) is ~ 1.3%. The largest fraction of dynamical mass 
of the host galaxy accreted through one satellite is 6.9% 
(10.1% for MRD profile) for G2 in NGC4526 (assum¬ 
ing that each GC structure corresponds to one accreted 
satellite). The largest fraction of total mass gained by the 
host galaxy through all the accretion events associated to 
the observed GC structures is ~ 12.2% (^ 17.8% for MRD 
profile) for NGC4526 (G1 and G2)). Among the other 
galaxies, the total fraction of dynamical mass contributed 
by all accretion events exceeds 4% only for NGC4365 
for the observed density profile; two more galaxies ex¬ 
ceed this limit when using the MRD profiles (NGC4365 
and NGC4621). The estimated dynamical masses of 
the accreted satellite galaxies span the 8.5, ~ 10.2] 
log M/Mq for the observed radial density profiles and 
[~ 8.3, ~ 10.3] log M/M© for MRD profiles. Considering 
that these are lower limits, they are consistently larger 
that the typical masses of dwarf elliptical galaxies (dEs) 


M ~ lO^M© ( 

Mateo et al. 

1993 

). Using the relation 

NgCs = obtainec 

[ from Harris et al. 

(2013) 


for dEs, the estimated masses of the accreted galaxies, 
on average, are only ~ 20% smaller than the masses ob- 
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Fig. 13.— Boundaries of the “large” (more than 19 pixels) GC over-density structures in the galactocentric distance (expressed in units 
of the effective radii re) vs azimuthal angle plane. The plots in the margins show the “large” structures in the R.A. vs Dec. plane. The 
background regions are color-coded based on their membership to the galaxies of our sample. 


tained with the relation for all elliptical galaxies, and 
still significantly more massive than typical dEs. This 
fact may represent another indication that these galaxies 
have undergone multiple mergers, which is a reasonable 
hypothesis since they are members of the Virgo clusters, 
and som e of them even in sub-groups in Virgo (see Sec¬ 
tion 


5.4). 


Figure compares the estimates of the distribution 
of dynamical masses of the satellites of the Virgo clus¬ 
ter galaxies (from Table 5), with the distributions of 
mass es of the Milky Way ( McMillan 2011) and M31 satel¬ 
lites (McConnachie 2012 H While we have used for all 
the M31 satellites the estimates of the dynamical m asses 


within the obser ved half-light radius available in (Mc- 
Connachle]|2012 1, different definitions have been used to 


determine the masses of the MW satellites, that usually 
yield mass estimates than can be systematically larger 
than the dynamical masses for the same systems. Even 
within the uncertainties and the possible systematic dif¬ 
ferences in the masses of the host galaxies, our estimates 
point to larger satellites in the Virgo galaxies. Consid¬ 
ering the similar mass of the MW, M31 and the Virgo 
ellipticals, environment may play a role. As shown, the 
mass of each accreted objects is a small fraction of the 
total mass of the main galaxy. Therefore, one would ex¬ 
pect only mild morphological and kinematic disturbances 
(assuming a dry merger) that would not be easily dis¬ 
cernible from the light of the galaxy. Furthermore, one 


would expect a correlation between the mass of the ac¬ 
creted satellite and the number of excess GCs. However, 
such a correlation is easily diluted by multiple cannibal¬ 
ism events, the time since coalescence, and the geometry 
of the interaction. Simulations of the evolution of the 
GC populations for a period of several Gyr after a dry 
merger are need ed to address these issues. 

In Section [5^ we reported that the GG structures ob¬ 
served in host galaxies located in lower-density regions of 
the Virgo cluster are more significant in the distribution 
of red GGs than in all GGs, than structures observed 
in galaxies located in higher-density regions. We argued 
that this evidence, in agreement with the predictions of 
cosmological simulation of the assembly history of galax¬ 
ies in cluster, indicates that recent major mergers are 
mostly responsible for the properties of the GG spatial 
distribution in the outskirts of Virgo, while the innermost 
galaxies, recently, have only accreted smaller companion 
galaxies. The fraction of dynamical mass contributed to 
the host galaxy by a single progenitor (shown in Table 5) 
does not exceed 5% (10% for MRD profiles) for the five 
galaxies (NGG4472, NGC4649, NGC4382, NGG4365 and 
NGC4621) located outside of the 33% isodensity contour 
shown in Figure [Tij No single GG structure accounts for 
25% of the mass of the final galaxy, which is the thresh¬ 
old used to label major mergers (i.e., mergers where the 
progenitors have 3:1 or lower mass ratio). Even assum¬ 
ing that all observed GG structures originate from a sin- 
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Fig. 14.— Positions of the ten Virgo cluster galaxies whose GC spatial distribution has been investigated in this paper (the D 25 isophotes 
of the galaxies, scaled for v isibility, are also plot ted). On the background, the surface density of confirmed and candidate Virgo cluster 
galaxies calculated from the |Binggeli et aL]||198^ catalog is shown. The two green lines are associated, from the center to the outskirts of 
the cluster, to the 50% (~ 40 galaxies/deg^) and 33% (~ 27 galaxies/deg^) fractions of the peak density. The insets display the location 
of the large GC structures within each host galaxy. 


gle accretion event, the total fraction of host galaxy dy¬ 
namical mass contributed by the progenitors (including 
also GC structures marked with asterisks in Table 5) 
for these five galaxies are 1.7% for NGC4472, 2.2% for 
NGC4649, 3.7% for NGC4382, 5.9% for NGC4365 and 
4.6% for NGC4621 (with MRD profiles the fractions are 
1.2%, 19.2%, 1.6%, 6.5% and 7.9% respectivel y). The 
disagreement between our discussion in Section |5.4| and 
these estimates depends mostly on two factor. Firstly, 
in our simple model that employs either the observed 
or MRD profiles to calculate the excess number of GCs 
for each GC structure, the dynamical evolution of the 
GCs is not taken into account. The changes in the mass 
function of GCs due to internal mechanism (stellar evo¬ 
lution and two-body relaxation) and to interactions with 
the host galaxy potential (gravitation al shocks and dy¬ 
namical friction; Fall & Zhang] 20011 all contribute to 
low er the estimated num ber ot excess GCs and, based on 
the Harris et al. (|2013 ) relations, to underestimate the 
dynamical mass oi the progenitors. The second factor is 
represented by the fact that major and/or wet mergers 
are usually conducive of new star formation and, con¬ 
sequently, the form ation of new GCs. The use of the 
relations derived by Harris et al. (20131 implies a “pas¬ 
sive” scenario where GCs are doiiated from the progeni¬ 
tors to the main galaxy; it does not take into account the 


effect of the merger on the GC system of the resulting 
galaxy. Moreover, the uncertainties on the GC specific 
frequencies for galaxies of large luminosity and differ¬ 
ent morphology that are typically involved in wet/major 
mergers and that can affect the scatter an d accuracy of 
the N nn - relation (see discussion inlHarris et al.f 

20131, further undermines the validity of the 


e estimated 

dynamical masses in the massive tail of the mass distribu¬ 
tion of the progenitors. It is also worth mentioning that 
major mergers are not the only possible way of forming 
new GCs, as wet mergers with mass ratios larger than 3:1 
may also produce significant amounts of star formation 
and, in turn, new young, metal-rich GCs. 

Detailed simulations of the formation and evolution of 
the GC populations in mergers of different types would 
be crucial to reconcile the tensi on b etween the assem¬ 
bly histories discussed in Section |5.4| and the inferences 
about the progenitor galaxies made in this Section, and 
would contribute to paint a more realistic and complex 
views of the evolution of GCs and their host galaxies in 
environment like the Virgo cluster. 

5.6. Caveats 

Our inferences on the pr ope rties of the progenitors of 
the GC structures (Section 5.5) are based on the hypoth¬ 


esis that all the structures originated by either accretion 
























































































TABLE 5 

Properties of the accreted galaxies associated to the “medium” and “large” GC over-density structures. 
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Fig. 15.— Comparison of the distribution of dynamical masses of the satellites of the Virgo galaxies estimated with the method disc ussed 
in this paper, with the distributions of masses of the satellite galaxies of M31 and the Milky Way collected by |McConnachie| (|2012f. The 
two boxplots on the right represent the distribution of inferred masses of the Virgo satellite galaxies assuming tnat tne UU over-density 
structures are due to accretion events. These are determined with the two relations and = ^0^^7±0.09 determined 

by|Harris et al.|(|2013|l using all elliptical galaxies of the sample and only dEs, respe ctively. The gray areas ab ove the boxplots represent the 
ranee ot mass valu es for the host ga laxies (the masses estimates for M31 are fromlEvans &c Wilkinsonl 2000l the interval of masses for the 
Milky Way is from [McMillan] | |2011[ l, and the single Virgo cluster galaxies are dynamical masses trom jttarns et al. |2013| l . The estimated 
masses of the individual satellites tor the Virgo cluster galaxies are represented as points with the same symbols and color-coding used in 
the other plots of th e pa per. Arrows are used for GC over-density structures that are likely to have been produced by major mergers (as 
discussed in Section |5.1[ |, and represent lower limits of the masses of the galaxies involved. The horizontal displacement of the points is 
only shown for the sake of clarity. 
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of satellite galaxies or major mergers. While this may 
be is an oversimplification of the multiple, diverse mech¬ 
anisms o f GC formation at work in massive early-type 
galaxies (|Brodie & Strader||2006|), it may apply at least 
for s ome of the (d(J structures we have detected (see Sec¬ 


tion 5.2). Below we discuss the possible issues that may 
affect our conclusions. 

5.6.1. Projection ejjects 

A concern is the possibility that the observed GC struc¬ 
tures are not physical structures, but may results from 
superposition of GCs, due to projection effects. Such 
bias cannot be excluded a priori without knowledge of 
the real three dimensional distribution of GCs around 
the host galaxy. It is nonetheless reasonable to assume 
that the importance of projection effects incre ases to¬ 
wards the center o f the galaxy. For instance, |Wood-| 
ley & Harris (2011) investigated the presence of phys¬ 
ical groups of GCs in NGC5128, a giant local elliptical 
galaxy, and argued that the chance of projection yield¬ 
ing chance two-dimensional groupings of all GCs is high 
within 10' (^2re). In our sample, because of the limited 
ACSVCS spatial coverage (see Figure [T^, we cannot ex¬ 
clude that some of the GC structures may be caused or 
enhanced by projection effects. For instance, streamer- 
like structure can be generated by three-dimensional GC 
spatial distributions following the triaxial stellar mor¬ 
phology shown by some elliptical galaxies. However, the 
risk of detecting unphysical GC structures is mitigated 
by the fact that most of the large and medium GC struc¬ 
tures at distance d < 2re have complex morphologies that 
cannot be explained with simple projection effects, as¬ 
suming a smooth three-dimensional spatial distribution 
of GCs. Moreover, pure projection effects cannot explain 
the striking differences observed in the properties of the 
same GC structures detected in distinct color and lumi¬ 
nosity classes (see Section |3 , even when we take into 
account the differences in the radial profile of red and 
blue GCs. 


5.6.2. Dynamical times 

Another possible issue with our interpretation of the 
origin of GC structures observed is related to their dy¬ 
namical evolution. The GC structures located in the 
outer regions of their host galaxies (A1 and A4 in 
NGC4472, Cl in NGC4649 and HI in NGC4365) might 
be relatively long-lived, as relaxation times are longer at 
large galactocentric radii. However, several observations 
have uncovered evidence of currently ongoing or recent 
accretio n of satellites in the halos of these galaxies (see 
Section 5.2). The GC structures located closer to the 
center of the host galaxies, assuming that they are real 
physical over-densities, ought to be the results of recent 
events in the assembly history of the host galaxy in order 
to be still visible today, given that relaxation times are 
much shorter at small galactocentric radii. This simple 
reasoning rules out the possibility that these structures 
have formed in the outskirts of the galaxies and drifted 
towards their observed positions under the sole effect of 
the dynamical friction of the host galaxy stellar field. 

We can estimate the time required for a single GC or¬ 
biting the host galaxy to reach the observed (projected) 
galactocentric distance from an arbitrary distance un¬ 
der the only effect of dynamical friction. We assume a 


standard mass for a single GC Mac = IO^Mq, consis¬ 
tent with the observed a verage absolute magnitudes of 
GCs (Strader et al. |2012 ), relative stellar circular veloc¬ 
ity fcirc~100 kms“^, which is the typical measured val¬ 
ues for Virgo cluster galaxies at the largest radius w here 
reliable estimates are available (Cortes et al.||2015l and, 
as initial distance rj = Ire. The angular distances of the 
closest dwarfs were convert ed to physical distanc es us¬ 
ing the galaxy distances by|Blakeslee et al.| (|2009|). The 
dynamical friction time AffHc of the GC infall has been 
calculated using equation (7-26) of Binney & Tremaine 
(2008), which express the time required for a GC to reach 
the center of the host galaxy as: 


^fric — 


1.17 rfucirc 
In A GM 


( 3 ) 


The hnal distance rj of the GC is set to the observed 
galactocentric distance of the GC over-density, which is 
a lower limit of the real three-dimensional distance from 
the center of the galaxy. The Atfric times range in the 6- 
30 Gyr interval for all GC structures in our sample (the 
structures Al and A4 in NGC4472, Cl in NGC44649, 
HI, H2, H3, H4 and H5 in NGC4365 and LI and L2 in 
NGC4552 were excluded because located at galactocen¬ 
tric distances larger than 1 re). In the case of spatially 
extended structures, the final position has been chosen 
to coincide with the center of mass of the structures. The 
uncertainties a{Atfric) can be as large as ~ 50% of the es¬ 
timated values, with average error (T(Atfric) 22%. The 
large times required for these structures to have reached 
the observed positions under the effect of dynamical fric¬ 
tion alone are incompatible with the typically short re¬ 
laxation times in the center of the galaxy that would 
have disrupted the GC structures and prevented the ob¬ 
servation of the GC overdensity structures. Assuming a 
fixed deprojected galactocentric distance r = 0.8 Ve, the 
dynamical friction times are, on average, ~ 13% shorter 
(up to ~ 40% in one case), making at least some of the 
observed GC structures compatible with the effect of in¬ 
ternal dynamical evolution only. However, while drift¬ 
ing towards the center of the galaxy under the effect of 
dynamical friction, the orbits of the GC belonging to a 
structure would lose coherence due to the differential ef¬ 


fect of tidal forces. For instance, Prieto & Gnedin (2008) 
found that GCs formed in a satellite halo moving on an 
initial circular orbit within the satellite halo, deviate sig¬ 
nificantly from the original orbit in a ^ 10 Gyr timeframe 
(shorter decoherence times are expected if the satellite 
galaxy has been tidally disrupted before 5 Gyr). GCs 
formed in the main halo and scattered by accreted satel¬ 
lite halos would deviate signihc antly from their original 
circular orbits within ~ 6 Gyr (Prieto & Gnedin 2008) 
from their formation. For these reasons, internal dynam- 
ical evolution alone can be mostly ruled out as explana¬ 
tion of the presence of the GC structures near the center 
of the galaxies. 

A possible formation mechanism for these GC struc¬ 
tures is that they are the coherent remains of the GC 
systems of satellite galaxies falling towards the center of 
the main galaxy at a large relative velocity. Another 
possibility is that recent, local formation o f new GCs 
has oc curr ed as a result of a m erger. While |C6te et 
(|1998| andlPipino et al.|(]2007|) found that the radial dis- 
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tribution, specific frequency and metallicity of the GC 
systems of massive elliptical galaxies can be understood 
by requiring a large fraction of GCs to have formed in 
satellite galaxies and stripped by the accreting galaxy, 
the formation of young, metal-rich GCs during merg¬ 
ers is believed to be another signihcant factor s haping 


of the observed GC systems of elliptical galaxies ( Bekki 


et al. 2003 fBrodie & Strader 2006 Wang et al. 2013 T 
A detailed modeling of the dynamics of the GG struc¬ 
tures progenitors, which is out of the scope of this paper, 
should also take into account the physical mechanisms 
acting on the GCs during their crossing of the accreting 
galaxy potential. These mechanisms would introduce dif¬ 
ferential effects in the spatial distribution of the observed 
GCs depending on the relative geometry of the orbit fol¬ 
lowed by the GC system and the potential of the host 
galaxy. 


5.6.3. Phase-space counterparts of spatial GC structures 

The physical existence of over-density structures in the 
projected spatial distribution of GCs can be confirmed if 
they are also characterized by a specihc motion, inherited 
by the accreted host galaxy or the merger progenitors, 
which is distinct from the systemic bulk motion of the 


main host galaxy (e.g. IRomanowsky et al.| 

2012 

Coccato 

et al.|2013 |Blom et al.l56l2ap. 'I'he lack oi 

spectroscopic 


data and kinematic measurements that would confirm 
the nature of all the GC structures in our sample will 
only be overcome as new observations become available. 
However, we stress that the kinematic identihcation of 
GC structures is a sufficient but not necessary condi¬ 
tion to ascertain the existence of physical groups of GCs 
originally associated to an accreted satellite galaxy. The 
kinematic method can only select GC structures with 
motions signihcantly different from the systemic veloc¬ 
ity of the host galaxy along the line-of-sight direction. 
GC structures mostly moving transversally to the line-of- 
sight are missed. Our approach could detect structures 
no matter what the direction of motion. 


6 . SUMMARY AND CONCLUSIONS 

We have studied the 2D distributions of the GC sys¬ 
tems of the ten most massive Virgo galaxies, us ing the 
GC catalog s extracted from the ACSVCS data (Jordan 
et al. 20091, an d GC catalogs extrac ted fro m analogous 
AGb d ata by (|Strader et al.||2012|) and (|Blom et al. 
2012b). We have estimated the statistical signihcance 
of localized spatial structures in these GC 2D distribu¬ 
tions, usi ng the KNN method a ugmented by MC sim¬ 
ulations (D’Abrusco et al. 2013). We have applied this 
method to the 2D distribution all the GCs in each galaxy 
and also to the classes of red, blue, low-L and high-L GCs 
(Section]^. Our results can be summarized as follows: 


• GCs over-density structures are common in the 
sample of galaxies studied. We observed a total 
of 229 GC structures in the residual maps of the 
spatial distribution of all GCs in the ten galaxies, 
and 42 significant, spatially extended GC struc¬ 
tures classified as “medium” or “large” structures, 
i.e. structures formed by groups of 5 to 19 adja¬ 
cent pixels or formed by 20 or more pixels with 
total simificance larger than 4a, respectively (see 
Section^. Most of these GC structures are also 


visible in the residual maps for red/blue and/or 
high-L/low-L GC subsamples. The physical sizes 
of the large and intermediate GC structures range 
between ^ 1 to ~ 25 kpc. 

• The spatial distribution of the small over-density 
structures is radially and azimuthally homogeneous 
(between I'and 2'). However, with the present data 
we cannot exclude that some of these small struc¬ 
tures may be due to statistical noise. 

• The medium and large structures may be related 
to the accretion of smaller satellite galaxies. This 
effect is stronger when the area of the GC struc¬ 
tures is considered. However, the spatial coverage 
of our sample is far from complete at larger radii. 


• The large, coherent spatial structures, which tend 
to be aligned along the major axis of the host galax¬ 
ies and are clearly visible in the spatial distribution 
of red GCs, could be the result of major mergers , 
as predicted by simulations (see Bekki et al.|2002 ). 
The medium-size structures may be related to the 
accretion of smaller satellite galaxies. 


• The smaller spatial structures (accounting for only 
^ 14% of the total area of the GC over-density 
structures observed) are observed to be more likely 
to be located along the direction of the minor axis 
of the host galaxy than the large structures than on 
the major axis. No signihcant differences are ob¬ 
served in their color and luminosity distributions. 


Large, spatially complex GC structures, mostly 
aligned along the major axis of the host galaxies, 
are more likely to be found in low-density regions 
within the Virgo cluster than in the high-density 
areas. While the red large GCs structures are, 
on average, more signihcant than struct ures in the 
spatial distribution of all GCs (Section 5.1), this 
ehect is stronger in the lower-density environment, 
suggesting the relatively recent mergers have deter¬ 
mined the GC spatial distribution in these galax¬ 
ies. Additional observations will be necessary to 
rule out observational biases. 


Using the total number of excess GCs associated to 
the large and medium GCs structures as a proxy 
of the size of the whole GC system, we have calcu¬ 
lated the range of dynamical masses of the satel - 
lites using the relations from Harris et al. (2013). 
We found that masses of the satellites estimated 
with our method are signihcantly larger than the 
observed masses of the satellites in the Milky Way 
and M31, suggesting that the environment plays 
an important role in the assembly of galaxies by 
accretion of satellites. The estimated masses of the 
progenitors are too small to be compatible with 
assembly histories involving recent major mergers, 
especially i n th e case of the host galaxies discussed 
in Section |5.4| We argue that, in case of major 
mergers, these disagreements could be reconciled 
with a more detailed description of the evolution 
of GC systems and the relation between size of the 
GC system and the properties of the host galaxy. 
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We also note that non-equal mass wet mergers with 
larger mass ratios could also be responsible for sig¬ 
nificant GC formation. 


We have proposed a simple mechanism responsible for 
the formation of these GC structures that invokes ei¬ 
ther small accretion of satellite galaxies or, in the case of 
the largest, spatially extended structures, major merg¬ 
ers, and used it to derive fundament al pr operties of the 
potential progenitor galaxies (Section 5.5). 

Our results show that 2D over-density structures in the 
spatial distribution of GCs in massive early-type galaxies 
are a common phenomenon, but we are currently limited 
by the lack of deep, high spatial resolution data reach¬ 
ing the outskirts of the host galaxies. New HST obser¬ 
vations complementing the ACSVCS data in the outer 
regions of the Virgo galaxies are necessary to fully deter¬ 
mine the characteristics of the 2D spatial distribution of 
GCs. Moreover, GCs catalogs covering larger galactocen- 
tric baseline would be useful to compare our results with 
data from other tracers of the dynamics and evolution of 
the host galaxy, like GCs and Planetary Nebulae kine¬ 
matics. The events in the assembly history of the host 
galaxy that, according to our model, may have generated 
the observed GC structures, should also leave evident sig¬ 
natures in the dynamics of these systems. Unfortunately, 
kinematic data for homogeneous, large samples of GCs 
located within the inner region of the host galaxies that 
could validate our assumptions are difficult to obtain. 

Constraining the evolution of massive galaxies also 
requires synergy between observations and simulations. 


Fully cosmological simulations of galaxy mergers are now 
being extended to resolve GC-sized objects (Mayer 2014, 
private communication) and to determine the effect of 
hierarchical evolution on the spatial distribution of the 
GC systems. This new generation of simulations will 
provide a unique testbed for the formation and evolution 
of GCs with unprecedented resolution and physical de¬ 
tail, producing GC parameters directly comparable with 
the observations. The comparison of the features of the 
observed 2D GC structures with the simulated data will 
assess how the observed 2D GC structures can be used 
to infer past merging/accretion history of the host galax¬ 
ies. The study of the GC over-density structures and its 
comparison with cosmological simulations that include 
GC formation will represent a powerful, game-changing 
combination of tools to advance the understanding of 
how galaxies grow and evolve. 
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APPENDIX 

COMPARISON WITH OTHER DENSITY ESTIMATION METHODS 

In order to confirm the effectiveness of our method and its ability to detect over-density structures in the projected 
spatial distribution of GCs, we have compared the residual maps obtained with our KN N-based technique with residual 
maps derived with a multidimensional Ga ussian Mixture Model (GMM) technique (Dempster et al. 1977) and the 
Kernel Density Estimation (KDE) method ( Rosenblatt||1956 |Parzen|[l962 ). 


The GMM method approximates the density function of the observed points with a combination of multiple, indepen¬ 
dent, multi-dimensional Gaussian distributions, where ea ch component repr esents a different group or subpopulation. 


2012) of the GMM technique, written in R (|^ 


In this paper, w e have used the Mclust implementation (Fraley et al. 

Core team||2014 ). In this implementation, the number of mixture model Gaussian components can be either specified 
by the user or determined automatically by the algorithm with an EM approach. R is an open-source free statistical en¬ 
vironment developed under the GNU GPL (http://www.r-project.org). An application of the Mclust algorithm to 
the problem of the detection of substructures in the velocity distribution of a sample of galaxies is discussed by|Einasto| 

eFalldM^. 

The KDE method is a non-parametric density estimation technique which reconstructs the unknown density function 
of a random variable using a kernel with a typical smoothing parameter. This parameter, the bandwidth, is estimated 
locally and depends on the local topology of the distribution of points whose density function is reconstructed. 

We have applied the KNN, GMM and KDE method s for the determination of the density maps to the observed 
distribution of ACSVCS GCs in NGC4472 (see Section 4.1), and to a simulated distribution of GCs where distinct 
components are visible. Density maps of the two GCs distributions have been determined in both cases. The other 
steps for the calculation of the residual maps from the density of the observed GCs and the model simulated GC 
distribution were unchanged. 

Spatial distribution of GCs in NGC4472 

The results of the application of the three density estimation techniques to the spatial distribution of the observed 
GCs distribution in NGC4472 are shown in F igur e [T6l The upper left plot shows the residual maps obtained with 
the KNN method (AT = 9) described in Section 4.1 the upper right plot is the residual map obtained using the KDE, 
while the lower panels show the residual maps obtained using the Mclust method where the number of components 
of the mixture model is determined automatically by the algorithm (left), or fixed to the optimal valu^right). The 

® We have derived the residual maps for number of components 
varying between 3 and 30, and selected the map the most closely 
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Fig. 16. — Residual maps of the GC distribution in NGC4472 obtained with different density estimation methods: upper left) KNN with 
K = 9; upper right) KDE; lower left) Mclust with number of mixture model components automatically determined; lower right) Mclust 
with optimal number of components fixed to 11. 

four residual maps are in substantial agreement on the overall number and positions of the GC structures detected by 
the KNN method, while the spatial extension varies significantly. The KDE and Mclust methods produce “smoother” 
and more spatially extended structures than the KNN method. 

Simulated spatial distribution of GCs 

We have also applied the KNN, KDE and Mclust methods to a simulated distribution of GCs (Figure!^ consisting 
in three different components: an azimuthally smooth distribution of GCs with radial positions randomly drawn from 
a power-law peaking in the center of the galaxy; three clusters, modeled as two-dimensional Gaussians with different 
variances along the R.A. and Dec. axes, radial and azimuthal positions relative to the center of the galaxy, and 
one “streamer”, modeled as a uniform azimuthal and radial distribution of GCs within a circular sector with 40deg 
aperture. The total number of simulated GCs is 1750, split in 1500 GCs belonging to the smooth distribution, 100 
GCs in the “streamer” and 50 GCs for each cluster. The total number of GCs and the number of GCs in each 
spatial component have been set to yield average GC densities and GC density contrasts in the over-density structures 
comparable to the observed values in the Virgo galaxies investigated in this paper (see Table for comparison). As 
a result, the number of excess GCs per pixels in each spatial structure of the simulated distribution is similar to the 
average number of excess GCs in the structures observed in the galaxies investigated. The spatial distribution of the 
simulated components are color-coded in the upper left panel in Figure [l^ As in the case of the NGC4472 GC system, 
the residuals maps produced with the three methods are similar: the two clusters E and N-W of the center of the 
galaxy are recovered in each map as well as the streamer (which, in the case of the KDE, has a rounded shape reflecting 
the use of a Gaussian as kernel function). The third cluster located just S of the center of the galaxy is associated to 
lesser over-density structures in the KNN and Mclust residual maps but it is missed by the KDE. Both the KDE and 
Mclust maps show artifacts structures located in the corners of the field and due to border effects, which are missing 
in the map produced by the KNN. 


Conclusions 

The three methods recover the same GC structures in both examples, with similar spatial and statistical properties. 
The shapes of the over-densities reconstructed with the GMM and KDE techniques, in general, have smoother bound¬ 
aries than those of the KNN-determined structures because the two former methods tend to spread the small-scale 
over-densities (often driven by GC excesses in one or few pixels) over larger areas. Moreover, GMM and KNN “bridge 

resembles the KNN residual map 
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Fig. 17.— Scatterplot of a simulated distribution of GCs with three distinct components (upper left) and residual maps obtained with the 
K = 9 KNN method (upper right), the KDE (lower left) and the Mclust method with number of mixture model components automatically 
determined (lower right). 

over” sharp pixel-by-pixel changes in the GC density, and the reconstructed final shapes are not representative, on spa¬ 
tial scales of few pixels, of the real geometry of the st ruc ture s. T he significances of the structures have been measured 
independently in the residual maps shown in Figures [T6| and [TT] as described in Section Overall, the significances of 
all GC structures in both GC distributions differ by less than l.Scr and O.Scr respectivdy, with negligible overall bias 
in the their distribution. The overall similarity of the significances indicates that the KNN method provides robust 
estimates of the real statistical significances of the GC structures. The small fluctuations measured are likely due to 
the different spatial extensions of the GC structures detected in the different residual maps. For instance, the lack of 
> 3(7 pixels in the map of the NGC 4472 GCs spatial distribution in Figure [T6]produced with the unconstrained Mclust 
algorithm is balanced by the larger areas occupied by the structures. The combined effect of singe pixels significance 
and total extension and shape of the structures produced similar statistical significances for all methods tested. 

Our conclusions is that, overall, the properties (location, significance and shape) of the GC residual structures 
obtained with the KNN method, are robust and in substantial agreement with the features of the GC structures 
detected using different density estimation methods. 
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